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EXECUTIVE  SUMMARY 


Introduction;  Both  satellite  radio  cosununlcatlons,  which  travel 
through  the  Earth's  Ionosphere,  and  high  frequency  (HF)  sky  wave  radio 
circuits,  which  use  the  Ionosphere  as  a  refracting  medium,  can  be  strongly 
affected  by  radio  wave  scintillation.  By  strongly  disturbing  the 
Ionosphere  high  altitude  nuclear  explosions  (HANEs)  can  cause  scintillation 
(and  other  effects),  thus  severely  degrading  satellite  (and  HF)  radio 
communications  over  a  large  region.  Iherefore  an  accurate  assessment  of 
Che  severity  of  nuclear  communications  effects  Is  necessary  In  order  to 
deploy  a  communications  system  which  Is  adequate  for  projected  needs,  but 
not  over-designed  and  hence  more  expensive  than  necessary.  Since 
experimental  evidence  Is  limited  to  a  few  high  altitude  nuclear  tests 
(before  the  test  ban  treaty  of  1963),  such  an  accurate  assessment  must  rely 
heavily  on  computer  simulations  and  thus  on  a  thorough  tmderstandlng  of  the 
physics  Involved  in  both  the  disturbed  ionosphere  and  its  Interaction  with 
radio  communications  signals.  Fortunately,  the  propagation  of  radio  waves 
through  a  given  disturbed  Ionosphere  is  relatively  well  understood. 
Accordingly,  our  1980  JASON  Summer  Study  focussed  primarily  on  some  aspects 


Insufficient  to  answer  many  questions  which  have  arisen  subsequently* 
Since  no  further  high  altitude  testing  Is  permitted  under  the  treaty,  one 


must  simulate  the  effects  of  high  altitude  nuclear  explosions  using  a 

I 

rather  elaborate  sequence  of  cooq>uter  codes*  These  codes  (illustrated  in 
Fig*  2  of  the  main  report)  are  based  on  current  understanding  of  the 
physics  Involved  In  high  altitude  nuclear  effects*  However,  In  some  areas 
the  physics  is  not  well  understood  and  simulations  must  be  checked  by 
experiment*  Such  experiments  have  Included  rocket  releases  of  barium  In 
the  ionosphere  and  observations  of  natural  Ionospheric  irregularities*  In 


this  report  (section  IV)  we  discuss  how  scaled  down  laboratory  experiments 
could  be  used*  Comparison  with  experiment  Is  the  key  way  to  check  whether 


simulations  have  overlooked  any  physical  processes* 


Current  Issues  Regarding  the  Physics  of  Ionospheric  Nuclear 


Effects!  The  theoretical  and  empirical  methods  for  forecasting  nuclear 
effects  currently  Involve  a  number  of  critical  issues*  A  key  factor  in 


these  Issues  is  the  enormous  range  of  scale  sizes  involved  In  the 


problem*  He  are  concerned  with  size  scales  as  small  as  ~  10  to  100  m,  yet 
the  region  of  Interest  is  thousands  to  tens  of  thousands  of  km  in  size— a 
size  range  of  ~  10^  to  10^*  Computer  simulations  can  handle  size  ranges  of 
about  10  to  100!  Perhaps  the  most  important  of  these  issues  is  the  physics 
which  governs  the  evolution  of  ionospheric  plasma  irregularities  at  late 
times  (~  minutes  to  hours  after  a  nuclear  event). 


iv 


O 


Do  the  Irregularities  evolve  for  a  few  linear  growth  times  of  the  relevant 
instability  (~  tens  of  minutes  for  a  large  barium  cloud)  and  then  “freeze" 
or  do  they  continue  to  evolve?  What  are  the  roles  of  currents  flowing 
along  geomagnetic  field  lines  and  connecting  across  the  field  in  a  remote 
region,  of  Alfv6n  waves,  of  finite  Larmor  radius,  etc*?  Bow  can  coiq>uter 
simulations  cope  with  the  enormous  range  of  scale  sizes  Involved?  We 
suggest  a  “fractal"  model  to  help  in  coping  with  some  of  these  questions. 
With  regard  to  the  question  of  i-iitial  seed  strlations,  we  suggest  that 
laboratory  experiments  using  laser  exploded  pellets  in  large  vacuum  tanks 
might  be  helpful. 

Other  Issues  involve:  (1)  the  possible  existence  of  plasma  which 
was  unobserved  optically  during  high  altitude  nuclear  tests,  but  which 
would  still  disturb  radio  wave  propagation  over  a  volume  substantially 
larger  than  that  of  the  optically  observed  plasma,  (2)  the  possibility  of 
magnetospherlcally  trapped  ions  forming  an  artificial  ring  current  and  the 
litq>act  of  such  a  ring  current  on  striated  plasmas. 
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Conclusions,  Recommendations  and  Action  Items;  Our  conclusions 
and  recommendations  can  largely  be  embodied  in  the  following  list  of  action 
items  through  which  the  recommendations  can  be  implemented. 

1.  Investigate  the  fractal  model  for  plasma  cloud  evolution  by 
(a)  finding  out  whether  or  not  the  evolution  of  ionospheric  barium  releases 
(as  we  now  understand  them)  can  be  viewed  as  fractalizatlon  processes, 
i.e.,  compare  observations  with  the  implications  of  a  fractal  model  and 


(b)  investigating  whether  or  not  the  equations  now  used  to  describe  plasma 
Irregularity  evolution  lead  mathematically  to  a  fractal  model*  The  fractal 
model  is  important  because  it  can  extend  existing  computational  models  to 
cover  the  million-fold  size  range  necessary  to  study  nuclear  environment 
scintillation  at  late  times  and  because  preliminary  results  of  the  fractal 
model  indicate  that  late  time  scintillation  may  be  more  disruptive  to 
satellite  communications  than  currently  anticipated* 


2*  Do  short  design  studies  to  find  out  id^iether  or  not  meaningful 
laboratory  experiments  (involving  laser  exploded  pellets  in  large  vacuum 
tanks  containing  magnetic  fields)  can  be  done  and  for  what  cost*  Such 


experiments  might  provide  information  on  early  time  (<  minutes)  coupling  of 
the  expanding  debris  cloud  to  the  atmosphere,  the  influence  of  early  time 


structure  on  late  time  irregularity  evolution  and  inner  and  outer  scale 


sizes  of  late  time  plasma  irregularities* 


during  the  next  large  barium  cloud  release*  Because  of  ’^adiative  transfer 
effects  (Horak  and  Whittaker,  1980)  observations  of  large  barium  clouds  can 
be  misleading  if  only  the  strong  resonance  transition  at  A55A  A  is 
observed* 


A*  Theoretically  investigate  scintillation  effects  alone  a 


typical  HF  skywave  path  through  a  late  time  nuclear  environment*  If  HF 


systems  are  to  provide  a  viable  backup  communications  mode,  it  is  necessary 
to  assess  the  severity  of  scintillation  effects  upon  them*  Nuclear  effects 
computer  codes  for  HF  propagation  do  not  now  Include  scintillation  effects* 


vl 


5.  Investleate  the  likelihood  and  consequences  of  an  artificial 
energetic  Ion  (ring)  current  being  formed  as  a  result  of  a  high  altitude 
nuclear  explosion*  If  such  a  magnetospherlc  ring  current  were  to  be 
present  following  a  high  altitude  nuclear  explosion  (and  there  is  no  direct 
experimental  evidence  one  way  or  the  other),  there  could  be  substantial 
consequences  for  the  evolution  of  ionospheric  irregularities* 


6*  Incorporate  radio  and  radar  data  into  the  description  of 
nuclear  event  plasma  following  Fishbowl  (and  possibly  other)  high  altitude 


nuclear  tests,  particularly  regarding  spatial  distribution  of  the  plasma* 
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1.  INTRODUCTION,  CONCLUSIONS  AND  RECOMMENDATIONS 

; 

■Both  satellite  radio  communications,  which  travel  throvigh  the 
Earth's  ionosphere,  and  high  frequency  (HP)  sky  wave  circuits,  which  use 
the  ionosphere  as  a  refracting  medium,  can  be  strongly  affected  radio 
wave  scintillation.  High  altitude  nuclear  explosions  cause  scintillation 
(by  strongly  disturbing  the  ionosphere)  and  thus  severely  degrade  satellite 
radio  communications  over  a  large  region.  Since  further  atmospheric 
nuclear  tests  are  banned,  a  thorough  understanding  of  the  physics  Involved 
in  both  the  disturbed  ionosphere  and  its  Interaction  with  radio  waves  is 
necessary  in  order  to  design  radio  communications  systens  which  will 
operate  satisfactorily  in  a  nuclear  environment.  During  the  1980  JASON 
Summer  Study  we  addressed  some  aspects  of  the  evolution  of  ionospheric 
irregularities  following  a  high  altitude  nuclear  explosion — the  radio  wave 
propagation  theory  being  apparently  well  understood  for  the  satellite  link 
case.  In  particular,  we  have  worked  on  Irregularity  evolution  at  late 

Ct  V  \  '  .V  . 

times  hours)  after  an  explosion  and  the  impact  of  early  time 
irregularity  structure  on  late  time  evolution.  We  also  raise  the  question 
of  scintillation  effects  on  HP  sky  wave  communications 

In  this  first  section  of  our  report  we  put  our  Investigation  in 
perspective  by  first  discussing  the  physical  phenomena  connected  with  high 
altitude  nuclear  explosions  and  the  impact  of  these  phenomena  on  radio 
communications  systems  using  satellites.  We  then  discuss  the  methods 
currently  used  to  estimate  these  effects  in  the  absence  of  complete  data 
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from  atmospheric  nuclear  tests,  namely  computer  simulations.  There  are  a 

number  of  ways  some  aspects  of  these  simulations  are  or  can  be  checked  both  ^ 

experimentally  and  theoretically.  He  review  existing  methods  and  suggest 

some  new  ones.  In  checking  the  computer  simulations  a  number  of  critical 

Issues  have  arisen.  These  are  outlined.  Finally  we  sum  up  our  summer  ^ 

study  conclusions  and  recommendations. 

A.  Impact  of  High  Altitude  Nuclear  Explosions  on  Communications  ' 

via  the  Ionosphere 

1.  Phenomenology  of  High  Altitude  Nuclear  Explosions 

As  will  become  clear,  the  complicated  structure  of  the 
atmosphere  and  the  multitude  of  physical  processes  at  work  following  a  high 
altitude  nuclear  explosion  make  quantitative  description  of  the  nuclear 
environment  a  difficult  and  complex  task.  The  brief  description  which 

I 

follows  draws  on  Wittwer  {1980b  Ac),  Zlnn  et  al.  (1966)  and  Hess  (1968). 

As  a  beginning,  high  altitude  nuclear  phenomena  can  be  divided  Into  three 
time  regimes  according  to  the  principal  physical  processes  occurring,  to 
wit 

Fast  Regime:  0  <  t  <  few  seconds; 

HHD  Regime:  few  seconds  <  t  <  minutes; 

~  ~  ( 

Electrostatic  Regime:  minutes  <  t  <  hours. 

In  the  Fast  Regime  the  nuclear  device  explodes  with  debris  having 

(  I 

an  Initial  velocity  of  ~  1500  km  s""^.  The  kinetic  energy  In  this  debris  Is 
coupled  to  the  surrounding  air  by  X-ray  and  UV  radiation  which  Ionizes 
ambient  air  at  large  distances  (>  lOO's  of  km)  and  by  local  processes  such 


as  collisions  and  charge  exchange*  In  this  regime  the  burst  debris  expands 
to  about  one  cyclotron  radius  perpendicular  to  the  magnetic  field  (~  30  km 
for  the  Starfish  event)  and  much  larger  distances  along  the  geomagnetic 
field*  Ass3nnetrle8  In  the  device  disassembly  and  In  the  ambient  atmosphere 
aided  by  various  kinetic  plasma  Instabilities  tend  to  produce  structure  in 
the  expanding  debris  cloud  and  ambient  air  to  which  It  couples* 

Within  a  few  seconds  a  magneto-hydrodynamic  (MHD)  shock  wave  forms 
and  begins  to  expand  Initiating  the  MHD  Regime*  As  the  shock  'snow  plows' 
through  the  now  partially  Ionized  ambient  medium  it  loses  energy  locally  by 
collisions  and  charge  exchange*  In  the  downward  direction  this  'magnetic 
bubble'  Is  stopped  by  local  energy  losses*  If  the  burst  Is  at  a  high 
enough  altitude  (such  as  the  400  km  Starfish  event),  the  upward  notion 
cannot  be  stopped  by  the  very  diffuse  upper  atmosphere*  The  upward 
traveling  plasma  moving  through  the  geomagnetic  field  generates  an  electric 
field  E  ■  -  v  X  f  which  polarizes  the  cloud*  Polarization  and  conduction 
currents  flow  In  the  lover  portion  of  the  magnetic  bubble*  These  currents 
result  In  a  viscous-type  drag  oil  the  upward  motion,  eventually  halting  It 
when  the  bubble  Is  ~  500  km  In  radius  (for  Starfish)*  The  magnetic  bubble 
collapses  having  dissipated  a  large  fraction  of  the  debris  energy  over  a 
large  region  of  space  (~  500  km  radius  perpendicular  to  i  and  from 
southern  to  northern  conjugate  zones  along  the  geomagnetic  field)* 

Following  the  bubble  collapse  the  heated  atmosphere  begins  to  ejqiand, 
'heaving'  large  amounts  of  plasma  upwards  along  geomagnetic  field  lines 
into  the  magnetosphere  (see  Fig*  1).  This  process  continues  In  the  next 
regime*  In  addition  to  the  Irregularities  existing  from  the  fast  regime. 
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Figurt  1  Electron  density  contours  |cm“'*) 
lerge  amounts  of  plasma  stored  ii 
undisturbed  ionosphere  normally 


Rayleigh-Taylor  and  ^  field  line  curvature  Instabilities  Introduce 
further  structure. 

At  "late  times"  of  from  minutes  to  many  hours  large  amounts  of 
plasma  continue  to  be  stored  In  the  magnetosphere  as  Illustrated  In  Fig.  1 
for  t  ~  0.5  hr.  This  plasma  subsequently  descends,  continually  Injecting 
plasma  Into  the  Ionosphere.  This  curtain  of  plasma,  only  several  degrees 
wide  In  the  geomagnetic  east-west  direction,  extends  to  larger  distances 
north  and  south  along  the  geomagnetic  field.  The  principal  plasma 
processes  of  Interest  at  this  time  cause  this  period  to  be  called  the 
Electrostatic  Regime.  Here  the  geomagnetic  ?  field  Is  taken  to  be 
constant  so  that  current  systems  and  electric  fields  are  the  primary  agents 
affecting  plasma  transport.  As  shown  In  Fig.  1,  there  are  large  gradients 
In  plasma  density  perpendicular  to  the  geomagnetic  field.  These  gradients 
In  conjunction  with  strong  neutral  winds  (either  induced  by  the  explosion 
or  at  late  times  of  natural  origin)  cause  the  plasma  to  be  imstable  to  the 
i  X  i  gradient  drift  Instability  and  small  scale  10  m  to  lOO's  of  km) 
plasma  density  Irregularities  develop  on  a  time  scale  of  >  10  minutes, 
lasting  for  hours.  The  fact  that  these  irregularities  exist  over  such  an 
enormous  range  of  size  scales  (~  10^  to  10^)  makes  computational  simulation 
Inherently  difficult.  This  difficulty  Is  discussed  In  Section  III  below. 

2.  Propagation  Effects  on  Radio  Communications 
The  large  extent,  high  density  and  small-scale 
Irregularity  structure  of  the  plasma  curtain  shown  in  Fig.  1  make  it  a 
major  barrier  to  radio  communication  links  passing  through  the  curtain  from 
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ground  to  satellite  terminals-  HP  sky  wave  communications  links  passing 
through  the  disturbed  Ionosphere  are  also  seriously  affected*  For  the 
first  few  minutes  to  tens  of  minutes  after  the  nuclear  explosion  (depending 
on  geometry,  link  frequency,  and  so  forth)  plasma  densities  and  collision 
frequencies  are  high  enough  In  the  curtain  that  radio  waves  are  highly 
absorbed*  Thus  satellite  communications  links  passing  through  the  plasma 
curtain  are  Interrupted  by  excessive  signal  loss*  This  situation  Is  often 
referred  to  as  "absorption  outage." 


As  plasma  densities  decrease  and  the  absorption  outage  period 
ends,  the  major  radio  communications  effects  are  due  to  small  scale  size 
Irregularities  In  plasma  density  and  hence  In  radio  refractive  Index*  As 
discussed  briefly  above  and  In  more  detail  In  Sections  II  and  III  below, 
the  plasma  shown  In  Fig*  1  Is  unstable  and  tends  to  form  Irregularities 
extended  along  the  geomagnetic  field  and  having  scale  lengths  lO's  of  m 
to  lOO's  of  km  perpendicular  to  the  geomagnetic  field*  Because  these 
Irregularities  are  Intense  (variance  >•  mean)  along  the  ray  paths  connecting 
ground  and  satellite  terminals  and  have  scale  lengths  less  than  the 
appropriate  Fresnel  zone  size  (<  few  kilometers),  they  can  cause  Intense 
signal  amplitude  scintillation  lasting  for  hours.  For  phase  scintillation 
and  strong  scatter  amplitude  scintillation  scale  sizes  larger  than  the 
Fresnel  zone  are  also  Important*  Signal  scintillation  Is  a  general  term 
referring  to  some  or  all  of  the  following  more  specific  effects; 


•  Deep  fades  In  signal  power  (lO's  of  dB) 

•  Signal  phase  fluctuations  (lO's  to  lOO's  of  radians) 
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•  Frequency  selective  effects 

•  Fluctuations  in  pulse  propagation  (group  delay)  time 

•  Fluctuations  in  angle  of  arrival  at  the  receiving  antenna 
a  Decorrelation  of  satellite  signals  arriving  at  spatially 

separated  earth  terminals 

The  aforementioned  irregularities  in  plasma  density  also  affect  HF 
sky  vave  communications  which  propagate  at  dekameter  wavelengths  between 
earth-based  terminals  via  refraction  in  the  ionosphere.  The  effects  on 
communications  signals  are  similar  to  those  for  translonospheric  links 
mentioned  above,  Including  both  absorption  outage  and  signal 
scintillations.  However,  HF  links  generally  use  small  bandwidths  (~  20 
kHz)  and  must  deal  with  ionospheric  variability  under  normal  conditions. 

So  the  HF  effects  of  high  altitude  nuclear  explosions  may  in  a  relative 
sense  be  less  severe  than  effects  on  translonospheric  microwave  links. 


While  the  communications  effects  of  high  altitude  nuclear 
explosions  are  Indeed  severe,  they  can  be  circumvented  entirely  or  at  least 
mitigated  by  appropriate  system  design.  Some  straight-forward  techniques 
Include  the  following: 


V 

l> 

Increased  transmitter  power  and/or  antenna  gain 
Increased  operating  frequency 

Receiver  (MODEM)  design,  e.g.,  coding  schemes,  phase  tracking 
loops,  etc. 

Alternate  link  paths  to  go  around  the  disturbed  region 
TWO  or  more  antennas  spatially  separated 
Operational  procedures,  e.g.,  use  of  smaller  ban^ldths 
during  disturbed  periods 
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These  techniques  can  certainly  be  effective*  However,  they  significantly 
increase  the  cost  of  communications  systems*  Hence  an  accurate  assessment 
of  the  severity  of  nuclear  communications  effects  is  necessary  in  order  to 
deploy  a  communications  system  which  is  adequate  for  projected  needs,  but 
not  over-designed  and  hence  more  expensive  than  necessary.  Such  an 
accurate  assessment  requires  a  thorough  understanding  of  the  evolution  of 
plasma  irregularities  following  a  high  altitude  nuclear  explosion*  Much 
research  has  been  done  to  understand  this  evolution,  including  the  research 
reported  here* 

B*  Simulation  of  High  Altitude  Huclear  Effects  Using  a  Sequence 

of  Computer  Calculations 

The  last  series  of  high  altitude  nuclear  tests  by  the  United 
States  occurred  during  the  Fishbowl  series  of  1962.  Further  atmospheric 
testing  was  banned  by  treaty  in  1963.  Although  much  data  was  collected 
during  the  Fishbowl  series,  this  data  set  is  insufficient  to  answer  many 
questions  which  have  arisen  since  then.  Since  no  further  testing  is 
permitted  under  the  treaty,  one  must  simulate  the  effects  of  high  altitude 
nuclear  explosions  using  a  rather  elaborate  sequence  of  computer  codes. 
These  codes  are  based  on  current  understanding  of  the  physics  involved  in 
high  altitude  nuclear  effects. 


In  Figs.  2a  and  2b  we  illustrate  a  possible  simulation  scheme* 

The  figure  takes  up  the  calculation  after  device  disassembly*  The  fast 
regime  (0  <  t  <  few  seconds)  is  handled  by  a  combination  of  six 
interrelated  computer  codos.  KLYSMA  handles  the  coupling  of  device  debris 
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Ffgure  2a  Ffow  chart  for  computer  cod«  u$«d  Jn  the  simufitlon  of  high 

attitude  nuclear  explosion  (HANE)  events.  The  sequence  is  con 
tinued  in  Figure  2b.  After  Ossakow  (1980). 


LATE  TIME  CODES  OUTPUT  ^ 
CONTOURS  OF  AVERAGE  N(x). 
V  lx).  E  lx),  T  lx.,  J  lx). 

ETC.  LARGE  SCALE  SIZES 


MAGIC 


BLOCKS  INSIDE  THE 
DASHED  BOX  REPRE¬ 
SENT  CONCEPTUAL 
FUNCTIONS  RATHER 
THAN  EXISTING 
COMPUTER  COOES 


LINEAR  STABILITY  ANALYSIS  L 
ON  THE  ABOVE  STATE. 

IF  STABLE  NO  IRREGULARITIES 
IF  UNSTABLE  PROCEED 


MICRO  AND  FLUID  CODES  jj 
THAT  FOLLOW  THE  EVOLUTION 
OF  THE  INSTABILITIES  TO 
SATURATION.  (NON-LINEAR) 
SMALL  SCALE  STRUCTURES 
50  KM) 


SPECTRUM  OF  IRREGULARITIES 
AS  A  FUNCTION  OF  SPACE 
^  (k,^) 


PROP 

PROPAGATION  CODES 

Figure  2b  Continuation  of  flow  chart  for  computer  simulation  of  h>^  altitude 
nuclear  explosion  (HANE)  events.  After  Ossakow  (1978). 
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energy  to  the  ambient  plasma,  heating  ions  and  electrons.  The  SCORPIO  and 
CYLUHAT  codes  use  coupling  region  parameters  generated  by  KLYSMA  to 
calculate  fluxes  of  energetic  ions  and  UV  and  X-ray  photons  emerging  from 
the  coupling  region.  The  CYLUHAT  code  handles  the  blast  wave  and  magnetic 
bubble  expansion  beginning  the  MHD  regime  (few  seconds  <  t  <  minutes).  The 
particle  and  photon  fluxes  are  deposited  over  a  large  region  by  the  PRODEP 
and  UVXRDEP  codes  for  particles  and  photons  respectively.  SWITCH  merges 
the  data  from  the  aforementioned  codes  into  an  input  to  three  dimensional 
MHD  codes  such  as  MRHYDE/PHOENIX,  which  calculate  the  atmospheric  motion  of 
coupled  ion  and  electron  fluids  as  large  amounts  of  plasma  are  heaved 
upwards  into  the  magnetosphere.  At  this  stage  (Fig.  2b)  an  electrostatic 
code  MAGIC  follows  the  subsequent  motion  of  the  plasma  produced  by  the 
nuclear  explosion  as  well  as  neutral  atmosphere  disturbances  at  lower 
altitudes.  We  are  now  in  the  electrostatic  regime  (minutes  <  t  <  hours) 
where  the  plasma  distribution  in  space  is  evolving  relatively  slowly.  Here 
the  time  evolution  of  plasma  density  N(x),  velocity  ^(x),  temperature  T(x) 
current  ^(x)  is  calculated  on  a  ~  50  km  grid.  Since  we  are  interested  In 
irregularities  on  scales  <  50  tan,  further  calculations  must  be  performed  as 
shown  in  Fig.  2b.  At  this  point  in  the  sequence  we  move  from  existing 
computer  codes  to  conceptual  functions  which  do  not  yet  exist  as  running 
codes  in  the  sense  indicated  in  Fig.  2a.  At  each  of  the  widely  spaced 
(~  50  km)  grid  points  calculated  by  MAGIC  a  linear  stability  analysis 
(LSA)  calculation  would  Indicate  whether  or  not  the  plasma  is  unstable  to 
the  E  x  ^  gradient  drift  instability.  If  so,  a  small-scale  size, 
electrostatic  fluid  model  (HSFC)  of  the  plasma  would  be  used  to  follow  the 


+  I 

yield  a  description  of  the  irregularity  structure  #(k)  on  scale  sizes 
of  ~  50  km  and  less  which  Is  the  Input  needed  by  the  radio  wave  propagation 

I 

codes.  These  propagation  codes  ultimately  calculate  the  statistical  | 

I 

characteristics  of  radio  signals  transmitted  through  the  disturbed 
Ionosphere  resulting  from  the  nuclear  explosion.  These  signal  statistics 
are  In  large  measure  the  "bottom  line"  of  the  simulation  so  far  as  late 
time  radio  communications  effects  are  concerned. 

At  the  present  time  the  conceptual  functions  described  above  and 
In  Fig.  2b  are  Implemented  by  an  abbreviated  set  of  calculations  along  the 
lines  Indicated  above  but  using  only  a  restricted  set  of  parameters  rather 
than  the  full  output  of  the  MAGIC  code.  For  example,  in  the  plane 
perpendicular  to  i  the  outer  scale  Is  set  at  10  km,  the  mean  square 
electron  density  fluctuation  <(An^)^>  at  <n^>/it  ,  the  spectral  Index  a 
for  plasma  turbulence  at  a  »=  2  and  so  forth.  Calculations  made  on  this 
restricted  set  of  parameters  are  then  used  as  Input  to  the  radio  wave 
propagation  codes.  Multiple-level-three-dlmensional  codes  are  presently 
under  development. 

There  are  two  principle  difficulties  In  such  an  involved  set  of 
computations.  First,  the  physics  on  which  the  computing  algorithms  are 
based  must  be  correct  enough  to  provide  useable  answers.  Second,  the  sheer 
size  of  the  problem  coupled  with  limited  storage  and  computer  time  means 
that  variables  of  Interest,  like  plasma  density,  cannot  possibly  be 

calculated  with  fine  resolution  in  time  and  space.  Current  computational  ( 


capabilities  limit  the  range  of  temporal  or  spacial  scales  covered  by  a 
given  computation  to  a  factor  of  from  about  10  to  100.  Thus  the  MAGIC  code 
which  handles  a  region  of  space  covering  thousands  of  km  must  use  a 
computational  grid  with  a  spacing  of  ~  50  km.  Yet  we  are  ultimately 
interested  in  irregularities  up  to  5000  times  smaller.  Thus  this  "dynamic 
range"  problem  necessitates  approximations  in  the  calculations!  scheme. 

The  discussion  above  indicates  how  the  problem  is  currently  being  handled; 
namely,  small  scale  sizes  are  left  out  of  the  computation  until  the  linear 
stability  analysis  stage  (LSA  in  Fig.  2b). 


C.  How  Can  the  Computer  Simulations  be  Checked? 


;  %  , 
*S>: 


The  con^utatlonal  simulation  described  above  covers  an 
enormous  range  of  scale  sizes  in  both  time  and  distance  as  well  as  being 
complex  physically.  Clearly  such  a  complicated  simulation  requires 
experimental  verification.  Although  some  features  of  the  computer 
simulation  are  verifiable  using  data  from  the  Fishbowl  test  series,  many 
questions  remain  and  no  more  high  altitude  nuclear  tests  are  planned 
because  of  the  test  ban  treaty.  Nevertheless  some  aspects  of  the  computer 
simulation  can  be  checked  by  comparison  with  experiments  involving  the  same 
physical  processes  which  are  active  in  the  nuclear  case.  Such  e^eriments 
'.zlude  releases  of  barium  in  the  ionosphere,  observations  of  natural 
ionospheric  irregularities  and  scaled  down  laboratory  experiments.  Further 
the  algorithms  used  in  the  simulation  codes  can  be  examined  to  find  the 
impact  of  any  physical  mechanisms  %ihich  are  not  currently  Included. 

PlasiBas  can  be  treated  by  both  fluid  and  particle  models.  Nhilc  the 
particle  models  are  necessarily  limited  due  to  relatively  small  numbers  of 
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particles,  comparisons  of  results  from  both  models  are  possible  and  could 
be  used  to  check  the  fluid  models  now  being  used.  Below  we  discuss  each  of 
these  methods  for  checking  the  computer  simulations  mentioned  above. 


1 .  Barium  Releases  in  the  Ionosphere 

In  the  electrostatic  regime  (minutes  <  t  <  hours) 
probably  the  most  pertinent  experiment  short  of  an  actual  test  is  the 
release  of  some  kilograms  to  tens  of  kilograms  of  barium  in  the  ionosphere 
such  as  the  Secede,  Stress  (McDaniel,  ed.,  1978)  and  Avefrla  (Simons  et  al. 
1980  and  McDonald  et  al.  1980)  experiments.  The  barium  released  is  ionized 
either  in  the  release  process  or  by  solar  UV  radiation.  These  Ba  ions  can 
be  observed  optically  in  resonance  emission  lines  (e.g.,  4554  A).  However, 
the  resonance  lines  can  be  optically  thick  and  observations  in  optically 
thin  fluorescence  lines  (e.g.,  5854,  6142  and  6497  A)  provide  necessary 
additional  information  for  mapping  the  distribution  of  a  Ba  ion  cloud  in 
some  situations  (Horak  and  Whitaker,  1980).  In  Fig.  3  we  see  the  unstable 
plasma  formed  into  striatlons  oriented  with  the  long  dimension  along  the 


geomagnetic  field.  Such  barium  clouds  have  been  observed  in  situ  by  rocket 


probes  and  remotely  by  optical  wavelength  photometry,  backscatter  radars 
and  forward  scatter  radio  transmission  through  the  cloud.  Because  high 
altitude  nuclear  events  and  barium  releases  are  both  subject  to  the  £  x  B 
gradient  drift  instability,  these  bariinn  release  experiments  provide  much 
experimental  data  which  is  relevant  to  the  evolution  of  ionospheric  plasmas 
resulting  from  high  altitude  nuclear  explosions.  Although  the  barium 
releases  are  not  completely  analogous  to  the  nuclear  case,  they  provide  an 
excellent  opportunity  to  exercise  and  check  the  computational  algorithms 
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Figure  3  Striated  barium  plasma  orientad  with  long  dimension  along  the  geomagnetic 
field.  This  photograph  of  the  mature  ion  doud  was  taken  about  40  minutes 
after  the  release  of  48  kg  of  Be  at  about  180  km  altitude.  After  Boquist  (1977). 
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used  for  the  electrostatic  regime.  So  far  barium  plasmas  have  been 
observed  for  up  to  about  three  hours  after  release.  Simulations  by 
computer  have  carried  out  for  times  up  to  about  10  minutes  after  release 
and  are  in  substantial  argument  vd.th  the  observations  over  this  time 
period.  However,  simulations  have  not  been  carried  out  over  longer  time 
periods  and  there  are  questions  (discussed  below)  as  to  plasma  behavior  at 
these  late  times  (~  hours). 

2.  Natural  Ionospheric  Irregularities 

Natural  ionospheric  irregularities  occur  predominately 
at  equatorial  and  auroral  latitudes.  As  with  the  barium  plasma  clouds,  the 
natural  ionospheric  plasma  in  these  regions  is  subject  to  (among  others) 
the  same  instability  thought  to  cause  plasma  irregularities  during  the 
electrostatic  regime  of  high  altitude  nuclear  events.  Thus  a  thorough 
understanding  of  these  natural  ionospheric  Irregularities  provides  another 
opportunity  to  check  the  physics  used  to  construct  nuclear  event  simulation 
algorithms.  Fejer  and  Kelly  (1980)  provide  an  in-<iepth  review  of  the 
current  understanding  of  ionospheric  irregularities.  Perusal  of  this 
review  will  reveal  that  many  of  the  people  active  in  the  simulation  of  high 
altitude  nuclear  events  have  also  contributed  to  the  understanding  of 
ionospheric  irregularities. 

3.  Laboratory  Experiments 

Experiments  with  laboratory  plasmas  have  not  so  for  been 
particularly  relevant  to  the  simulation  of  high  altitude  nuclear 
explosions.  However,  we  argue  in  Section  IV  below  that  laboratory 


experiments «  specifically  involving  laser  exploded  pellets,  can  probably 
provide  important  information  on  the  following: 

•  Coupling  of  the  expanding  debris  cloud  to  the 
ambient  ionosphere  in  the  fast  and  MHD  regimes 

•  Importance  of  "early  time"  phenomena,  particularly 
spatial  structure,  on  subsequent  evolution  of  the 
plasma 

•  Inner  and  outer  scale  sizes  of  plasma 
Irregularities  in  the  electrostatic  regime 

Well  conceived  laboratory  experiments  could  provide  a  valuable  complement 
to  the  barium  release  and  natural  ionospheric  irregularity  experiments 
discussed  above  by  their  ability  to  explore  early  time  phenomena  and  to 
explore  late  time  effects  in  a  highly  controlled  environment.  Laboratory 
experiments  begun  at  this  time  would  benefit  from  the  large  scientific 
effort  devoted  to  controlled  thermonuclear  fusion  and  space  plasma  related 
laboratory  experiments. 

A.  Impact  of  Physical  Mechanisms  not  Currently  Included  in 
Computational  Simulations 

Computational  simulations  currently  used  to  simulate 
high  altitude  nuclear  explosions  must  deal  with  an  almost  incredibly 
complex  series  of  events  involving  a  large  number  of  physical  processes. 

In  order  to  make  the  problem  tractable,  only  those  effects  thought  to  be 
dominant  are  included.  In  order  to  make  the  simulation  valid,  any  effects 
not  Included  in  the  simulation  should  be  considered  and  either  Included  or 
shown  to  be  Insignificant.  For  example,  in  the  electrostatic  regime  finite 
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Larmor  radius  and  some  three  dimensional  effects  are  not  currently 
Included.  High  altitude  nuclear  explosions  inject  large  numbers  of  high 
energy  particles  into  the  magnetosphere,  some  of  which  are  trapped  for  a 
significant  period  of  time.  The  interaction  of  these  disturbances  of  the 
magnetosphere  with  the  ionosphere  is  not  now  included  and  may  be  important 
at  late  times. 


5.  Comparison  of  Fluid  and  Particle  Simulations  of  Plasma 
Dynamics 


Electrostatic  fluid  simulations  of  a  plasma  are  valid 
over  a  large  range  of  conditions.  However,  fluid  simulations  are 
macroscopic  and  thus  do  not  include  particle  effects  such  as  the  finite 
Larmor  radii  of  ions  and  electrons  in  the  plasma.  Particle  simulations 
such  as  those  by  Goede,  Dawson  and  Humanlc  (1980)  and  Goede  (1980)  can 
simulate  situations  which  are  relevant  to  the  evolution  of  ionospheric 
irregularities.  While  these  particle  simulations  cannot  at  present  attack 
the  ionospheric  irregularity  problem  itself  (because  of  computer  speed  and 
memory  limitations)  they  can  help  check  the  electrostatic  codes  now  used  by 
Kesklnen,  McDonald  and  Ossakow  (1980),  McDonald,  Ossakow  and  Zalesak 
(1978),  Keskinen,  Ossakow  and  Chatervedi  (1979),  Ossakow,  Zalesak  and 
Zabusky  (1977).  This  cross-checking  process  would  involve  finding  a 
relevant  domain  of  physical  parameters  which  both  methods  could  handle  and 
comparing  the  results  of  electrostatic  and  particle  simulations  for  the 
non-linear  evolution  of  an  unstable  plasma.  It  should  also  be  kept  in  mind 
that  fully  electromagnetic  particle  simulatic  ns  are  available. 
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D.  Current  Issues  Regarding  the  Hiyalcs  of  Ionospheric  Plasmas 
Resulting  from  High  Altitude  Nuclear  Explosions 


Above  we  have  discussed  theoretical  and  empirical  methods  for 
checking  existing  computer  simulations  of  high  altitude  nuclear  effects* 
This  checking  process  currently  revolves  around  a  number  of  critical 
Issues*  Below  we  outline  a  number  of  these  issues  which  are  of  In^ortance 
In  the  research  reported  here*  Our  conclusions  regarding  these  questions 
and  recommendations  for  investigating  these  Issues  are  summarized  In 
section  E  below. 


1.  Evolution  of  Plasma  Irregularity  Structure  During  the 
Electrostatxw  Regime  fmlnutes  <  t  <  hours) 


Both  experimental  studies  using  barium  releases  In  the 
Ionosphere  (McDaniel,  1978)  and  theoretical  studies  using  numerical 
simulations  of  non-linear  plasma  dynamics  (Ossakow  et  al*  1977)  Indicate 
that  Ionospheric  plasmas  unstable  to  the  ^  x  ?  gradient  drift  Instability 
will  tend  to  form  "fingers"  along  the  E  x  drift  direction  (see  Fig*  18) 
and  that  these  fingers  will  bifurcate  forming  further  smaller  fingers* 

After  a  few  linear  growth  times  (~  tens  of  minutes  for  a  large  barium 
release)  the  statistical  properties  of  the  irregularities  can  be 
characterized  by  a  one"dimenslonal  spectral  index  a  of  about  2  and  an 
outer  scale  size  (L^)  on  the  order  of  the  cloud  size  (see  Section  ID* 

The  question  at  Issue  Is  what  happens  next?  Does  the  cloud  continue  to 
bifurcate  and  evolve  with  time  or  does  it  become  stable  and  "freeze,"  l*e*. 


retain  the  same  statistical  properties  (a  and  )  ,  while  moving  In 


space,  dispersing  and  recombining?  Put  another  way,  what  physics  governs 
the  irregularity  evolution  at  late  times  and  what  eventually  limits 
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irregularity  evolution?  What  are  the  roles  of  currents  flowing  along 
geomagnetic  field  lines  and  connecting  across  the  field  in  a  remote  region, 
of  Alfv4n  waves,  of  finite  Larmor  radius,  etc.?  As  explained  above  and  in 
detail  in  Section  11  a  and  are  crucial  parameters  in  the  design  of 

satellite  communications  systems  for  use  in  a  nuclear  environment. 


Experimentally  and  theoretically  the  ability  to  answer  these 
questions  is  now  limited.  Barium  clouds  have  only  been  observed  for  at 
most  a  few  hours  and  usually  less,  and  computer  simulations  are  limited  to 
a  few  generations  of  bifurcation  by  memory  size  and  computation  time. 
Further,  one  should  continue  to  critically  examine  the  extrapolation  of 
knowledge  gained  from  surrogate  experiments  with  barium  clouds  to  the 
nuclear  case. 


2.  Origin,  Character  and  Impact  of  Fast  and  MHD  Region 
Structure  on  Irregularity  Evolution  at  Later  Times 


Observations  during  high  altitude  nuclear  tests  indicate 
that  small  scale  (<  km's)  structure  exists  at  very  early  times.  Since  the 
MHD  regime  simulations  (see  Fig.  2)  only  yield  data  on  a  ~  50  km  grid, 
smaller  scale  (<  50  km)  structure  has  not  been  considered  prior  to  the 
linear  stability  analysis  (LSA  in  Fig.  2b).  The  question  is,  then,  what 
sort  of  small  scale  structure  exists  at  early  times,  if  any,  and  how  does 
it  affect  the  structure  at  later  times?  Current  computations  beginning  at 
LSA  in  Fig.  2b  assume  a  given  structure  on  scales  less  than  ~  50  km.  How 
close  to  the  truth  are  the  assumed  "seed”  irregularities?  How  would 
alternate  choices  of  "seed"  irregularities  on  Important  size  scales 
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(<  km  s)  affect  the-  stability  of  the  plasma  and  the  non-linear  growth  of 
irregularities?  As  a  first  step,  Pongratz  and  Simone  (1979)  and  Simons  et 
al.  (1980)  have  addressed  the  impact  of  kinetic  plasma  instabilities  on 
this  problem  both  experimentally  (Avefria)  and  theoretically.  Of  course, 
such  linear  analyses  must  be  supplemented  1^  further  non-linear  analysis. 
Ossakow  (1980)  plans  to  address  the  problem  by  Inserting  "seed" 
irregularities  (on  size  scales  <  50  km)  into  existing  simulation  codes 

such  as  those  reported  by  Ossakow  et  al.  (1977).  We  argue  that  laboratory 
experiments  discussed  in  E  below  and  in  Section  IV  could  be  helpful  here. 

3.  Possible  Magnetospheric  Effects 

It  is  well  known  that  high  altitude  nuclear  explosions 
Inject  enormous  numbers  of  energetic  ions  into  the  magnetosphere  (Hess, 
1968).  While  conventional  wisdom  indicates  that  these  ions  are  cooled  by 
collisions  and  fall  back  down  into  the  ionosphere,  it  is  possible  that 
pitch  angle  scattering  could  establish  energetic  (>  1  keV)  ions  in  high 
altitude,  virtually  collisionless  orbits.  If  this  should  happen,  a  partial 
ring  current  (analogous  to  the  natural  ring  current  at  L  ~  4)  would  be 
established  lasting  from  ~  hours  to  days. 

Since  ground  based  optical  observations  would  not  reveal  such  ions 
and  no  relevant  satellite  data  from  past  high  altitude  nuclear  tests  exists 
(to  the  author's  knowledge),  we  can  neither  confirm  nor  deny  the  existence 
of  such  an  artificial  ring  current.  Further,  there  is  apparently  no 
theoretical  work  bearing  on  the  problem. 


Here  we  simply  raise  the  question  of  the  possible  impact  of  such 
magnetospheric  effects  on  the  scintillating  environment  following  a  high 
altitude  nuclear  explosion.  For  example,  this  artificial  ring  current 
would  presumably  have  sharp  gradients  which  would  drive  field-aligned 
currents  which  could,  in  turn,  interact  with  conventional  ionospheric 
strlatlon  mechanisms  by,  for  example,  driving  a  current  convective 
instability.  Electric  fields  parallel  to  the  geomagnetic  field  might  also 
be  generated  providing  cross-field  electron  transport  which  would  alter 
conventional  striation  processes  by  allowing  ion  diffusion.  A  modest 
theoretical  effort  seems  in  order  to  assess  the  likelihood  of  such  an 
artificial  ring  current  being  formed  and,  if  so,  what  its  likely  impact  on 
striation  physics  would  be. 

A.  Geophysical  Distribution  of  Plasma  Following  a  High 

Altitude  Nuclear  Explosion 

Much  effort  has  been  expended  in  trying  to  understand 
the  photographic  data  collected  during  the  Fishbowl  series  of  high  altitude 
nuclear  tests.  Work  by  Kilb  (1980)  and  others  has  documented  the  large 
scale  distribution  of  plasma  including  the  exlstance  of  jets  moving  across 
field  lines  and  deposition  of  plasma  in  unexpected  regions,  e.g.,  far  south 
of  the  Southern  Conjugate  Region.  Most  of  this  work  has  centered  on 
photographic  data  from  a  number  of  sites.  However,  large  amounts  of  plasma 
can  remain  undetected  by  optical  means,  especially  at  late  times.  Chesnut 
(1980)  mentioned  several  ground  and  air  based  radio  and  radar  observations 
during  Fishbowl  which  apparently  sensed  plasma  unobserved  optically. 
Examination  of  these  pieces  of  radio  and  radar  data  could  substantially 
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change  our  picture  of  where  and  over  what  time  scale  enhanced  ionospheric 
plasma  densities  and/or  Irregularities  exist  following  a  high  altitude 
nuclear  explosion. 

E.  Conclusions  and  Recommendations 

Belcv  we  summarize  the  conclusions  and  reconmendatlons 
arising  from  the  1980  JASON  Summer  Study  on  the  topic  of  Irregularities  In 
Ionospheric  plasmas  and  their  Impact  on  radio  communications.  Both 
satellite  radio  communications  links,  which  travel  through  the  ionosphere, 
and  HF  sky  wave  circuits,  which  use  the  Ionosphere  as  a  refracting  medium, 
can  be  strongly  affected  through  radio  wave  scintillation  caused  by  plasmas 
generated  by  high  altitude  nuclear  explosions.  Since  further  atmospheric 
tests  are  banned,  a  thorough  understanding  of  the  physics  involved  In  both 
the  ionospheric  plasma  and  Its  interaction  with  radio  waves  is  necessary  in 
order  to  design  radio  communications  systems  which  will  operate 
satisfactorily  In  a  nuclear  environment.  As  the  radio  wave  propagation 
theory  Is  apparently  »7ell  in  hand  for  the  satellite  link  case,  we  have 
directed  our  efforts  toward  understanding  the  physics  of  Ionospheric  plasi>a 
irregularities  following  a  high  altitude  nuclear  explosion  and  In 
particular  the  evolution  of  irregularities  at  3nte  times  (~  hours)  after  an 
explosion.  We  also  raise  the  question  of  scintillation  effects  on  HF  sky 
wave  communications. 

1.  Large  Dynamic  Range  and  Fractalization 

The  scale-free  equations  used  (e.g.,  by  Perkins  et  al* 
1973  and  Ossakow  et  al.  1977)  to  study  the  evolution  of  plasma 
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irregularities  in  the  electrostatic  (late  time)  regime  suggest  that  the 

irregularity  evolution  can  be  described  by  a  fractallzatlon  process  V 

(fingers  have  fingers  have  fingers,  etc.)  Using  the  fractallzatlon  model 

one  can  estimate  the  vei^  important  spectral  index  a  of  a  fully  developed 

irregularity  cloud  (see  Sections  II  and  III)  using  only  Information  ■■ 

available  from  a  computer  simulation  of  the  first  couple  of  fractalizations 

(bifurcations).  Since  limitations  of  time  and  memory  size  virtually 

prohibit  computer  simulations  for  more  than  a  few  generations,  the  fractal 

model  could  be  extremely  useful  in  estimating  the  evolution  of  plasma 

irregularities  at  late  times  (•''  hours)  after  a  high  altitude  nuclear 

burst.  Because  fractallzatlon  leaves  large  scale  features  substantially 

unchanged  as  time  progresses,  it  is  possible  that  a  plasma  cloud  is 

"frozen"  in  visual  appearance  as  it  drifts  across  the  sky,  but  is 

nevertheless  evolving  according  to  a  fractal  model  on  such  small  size 

scales  that  no  changes  are  noticed  in  overall  appearance.  Probably  most 

important  the  fractal  model  Implies  that  plasma  irregularities  could 

eventually  evolve  such  that  the  spectral  index  of  the  irregularities  ot 

becomes  substantially  less  than  2  (~  1.5).  This  result  is  unexpectedly  low 

and,  if  accurate,  would  imply  more  severe  communications  disruptions 

following  a  high  altitude  nuclear  explosion  than  are  now  anticipated! 

Because  it  applies  to  time  regimes  otherwise  inaccessible  and  has 
preliminary  results  which  are  surprising,  we  recommend  that  the  fractal 
model  be  pursuer*  further  both  in  terms  of  the  model  itself  and  its 
consistency  with  existing  knowledge,  l.e.,  can  the  evolution  of  a  barium 
cloud  as  we  now  understand  it  be  viewed  as  a  fractallzatlon  process? 
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2.  Laboratory  Experiments 


< 


I 


We  suggest  that  laboratory  experiments,  specifically 
laser  exploded  pellet  experiments  in  large  vacuum  tanks  containing  magnetic 
fields,  could  be  useful  in  obtaining  a  more  complete  understanding  of  high 
altitude  nuclear  explosions.  After  preliminary  examination  it  appears  that 
such  experiments  can  provide  information  regarding  the  following: 

a.  The  'early  time'  (fast  and  MHD  regimes)  coupling  of 
the  expanding  debris  cloud  to  the  ambient 
atmosphere 

b.  the  importance  of  early  time  phenomena,  e.g. , 
spatial  structure,  in  determining  phenomena  at 
subsequent  times 

c.  Ihe  inner  and  outer  scale  sizes  of  plasma 
irregularities  in  the  late  time  (electrostatic 
regime) 

These  experiments  would  impact  current  understanding  both  from  the 
experiutents  themselves  and  the  improvements  in  theory  and  computer 
simulations  which  would  follow  from  efforts  to  understand  the  experimental 
data.  Our  discussion  does  uot  constitute  an  experimental  design;  instead 
it  provides  a  strong  indication  that  relevant  experiments  can  be  designed. 


We  recommend  support  of  short  term  (~  3  month)  design  studies  by 
organizations  having  the  necessary  pieces  of  major  equipment,  e.g.,  high 
power  lasers,  large  vacuum  tanks,  etc.  If  these  design  studies  Indicate 
that  meaningful  experiments  can  be  performed,  we  recommend  that  they  be 
undertaken  with  provision  for  participation  by  the  appropriate  research 
groups  engaged  in  the  simulation  of  high  altitude  nuclear  effects. 
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3.  Impact  of  Scattering  by  Small~Scale  Ionospheric 

Irregularities  on  High  Frequency  Sky  Wave  Communication 

Although  satellite,  cable  and  microwave  relay  links 
carry  the  bulk  of  U.S.  military  and  civil  communications,  point  to  point 
communications  by  HF(3-30  MH2)  sky  wave  (ionospheric  refraction) 
propagation  are  frequently  the  best  (and  often  the  only)  backup  systems 
available  for  long  distances.  If  HF  sky  wave  systems  are  to  provide  viable 
backup  communications  when  satellite  links  are  unavailable,  they  must 
function  when  the  ionosphere  is  disturbed  by  high  altitude  nuclear 
explosions*  While  large  scale  nuclear  effects  on  sky  wave  propagation  have 
been  investigated  extensively,  the  effects  of  signal  scintillation  by 
small-scale  irregularities  have,  to  the  author's  knowledge,  been  ignored* 
Since  small-scale  irregularities  will  prestimably  exist  for  hours  after  a 
high  altitude  nuclear  explosion  and  since  natural  irregularities  of  similar 
characteristics,  though  probably  less  Intense,  severely  disrupt  sky  wave 
communication  in  the  equatorial  and  auroral  regions,  we  recomm-'jnd  that  the 
impact  of  small-scale  irregularities  on  HF  skywave  propagation  should  be 
assessed.  If  this  Impact  is  significant,  which  we  anticipate  it  will  be, 
we  recommend  it  be  Included  in  the  computer  simulations  of  HF  propagation 
in  a  nuclear  environment*  ihlc  would  enable  communication  system  designers 
to  mitigate  these  scintillation  effects  by  appropriate  design 
modifications* 

While  small-scale  ionospheric  Irregularities  are  a  nuisance  so  far 
as  HF  communications  are  concerned,  they  make  possible  VHF  (~  30  to  50  MHz) 
ionospheric  scatter  propagation  over  1500  km  distances*  Ionospheric 
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scatter  communications  utilizing  natural  Irregularities  were.  In  fact,  used 
to  affect  communications  between  the  U.S,  and  Europe  (via  Canada, 

Greenland,  Iceland,  and  Britain)  prior  to  satellite  systems* 

4.  Radiative  Transfer  Effects  In  Optical  Observations  of 

Barium  Ion  Clouds 

The  barium  Ion  has  two  resonance  transitions  at  4934  and 
4554  A  with  three  fluorescent  transitions  at  slightly  longer  wavelengths 
5854,  6142  and  6497  A.  In  the  past  the  common  procedure  has  been  to 
observe  only  the  strong  resonance  line  emissions  at  4554  A  (and 
occasionally  the  other  member  of  the  pair  at  4934  A)*  We  have  studied  work 
still  In  preparation  by  Horak  and  Whitaker  (1980)  and  agree  with  their 
conclusion  that  observations  In  the  resonance  lines  alone  can  yield 
erroneous  values  of  the  column  density  of  barium  Ions  and  that  observations 
can  be  unambiguously  Interpreted  if  supplemented  by  simultaneous  photometry 
of  the  fluorescence  lines*  We  also  note  that  strontium  is  an  almost 
unavoidable  contaminant  of  barium  and  therefore  neutral  Sr  Is  likely  to  be 
present  In  Ionospheric  barium  clouds*  The  neutral  Sr  has  a  resonance  line 
at  4607  A  which  could  contaminate  the  4554  A  barium  Ion  line 
measurements*  This  again  argues  for  redundant  spectral  observations  (both 
resonance  and  fluorescence  lines)  in  cases  where  large  optical  depths  are 
anticipated — the  forthcoming  PLACES  experiment  should  produce  optical 
depths  of  ~  50  and  saturation  effects  In  the  resonance  line  at  4554  are 
certainly  likely*  We  therefore  recommend  that  both  resonance  and 
fluorescent  line  observations  be  made  during  PLACES  and  other  future  barium 
releases  where  large  optical  depths  at  4554  A  are  anticipated* 
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5.  Possible  Magnetospheric  Effects 


High  altitude  nuclear  explosions  inject  large  numbers  of 
energetic  Ions  Into  the  magnetosphere.  Some  of  these  Ions  may  become 
trapped  for  considerable  periods  (~  hours  to  days)  in  colllslonless 
orbits.  If  such  an  artificial  ring  current  were  to  be  present  following  a 
high  altitude  nuclear  explosion,  there  could  be  substantial  consequences 
for  the  evolution  of  ionospheric  irregularities.  A  modest  theoretical 
effort  seems  in  order  to  assess  the  likelihood  of  such  a  partial  ring 
current  being  formed  and,  if  so,  what  its  likely  impact  on  striation 
physics  would  be. 


i 


II.  IMPACT  OF  PLASMA  IRREGULARITIES  ON 
CCmiNICATIONS  VIA  THE  IONOSPHERE 

Since  the  refractive  index  for  electromagnetic  waves  varies  from 

unity  in  even  a  weakly  ionized  plasma,  the  large  amounts  of  ionospheric  and 

magnetospheric  plasma  generated  hy  high  altHude  nuclear  explosions  (Figure 

1)  are  likely  to  have  a  significant  influence  on  radio  communication 

signals  propagating  through  these  ionized  regions.  For  radio  signals  of 

interest  here  the  refractive  index  of  a  simple  ionized  medium  is  given  by 

n'*  •  1  -  (fp^/f  )  where  f  is  the  wave  frequency  and  the  electron  plasma 

frequency  is  f  «  10^  ^/N' where  N  is  the  number  density  of  electrons  in 
pe  ' 

cm'^.  Thus  one  finds  that  (other  factors  being  equal)  communication 
effects  are  more  pronounced  at  lower  radio  wave  frequencies  and  at  hi^er 
plasma  densities.  These  effects  apply  both  to  signals  propagating 
completely  through  the  nuclear  disturbed  ionosphere  and  magnetosphere, 
e.g.,  to  communications  satellites,  and  to  signals  propagating  via 
ionospheric  refraction  between  two  surface  locations,  e.g.,  sky  wave  HF 
communications.  At  times  not  long  after  a  high  altitude  nuclear  explosion, 
i.e.,  <  tens  of  minutes,  the  nuclear  generated  plasma  is  so  dense  and  so 

strongly  coupled  to  the  ambient  neutral  atmosphere  that  absorption  effects 
dominate  and  received  radio  signals  are  simply  too  weak,  too  greatly 
attenuated,  to  be  useful.  As  time  passes,  absorption  effects  decay  and 
scattering  effects  dominate,  lasting  hours.  In  assessing  scattering 
effects  the  spatial  size  distribution  of  the  plasma  irregularities  is  of 
great  importance.  The  evolution  in  time  of  the  spatial  distribution  of 
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nuclear  produced  plasma  Is  the  primary  subject  of  this  report.  Hence  we 

,1 

\ 

will  deal  with  propagation  and  communications  effects  only  long  enough  to 
indicate  the  motivation  for  studying  the  plasma  spatial  structure. 

( 

A.  Summary  of  Communications  Effects  Caused  by  Plasma 
Irregularities 

The  propagation  of  UHF  and  microwave  communications  signals 
(  >  300  MHz)  through  the  plasma  generated  by  a  high  altitude  nuclear 
explosion  has  been  studied  in  great  detail  by  Rlno,  Vondrak,  and  Hatfield 
(1975),  Johnson,  Chesnut,  and  Rlno  (1976),  Wlttwer  (1976),  Wittwer  et  al. 

(1977),  Hendrick  (1977),  Knepp  (1977),  Wittwer  (1978),  McDaniel  (1978)  and 
Wittwer  (1980a)  among  others.  Although  refractive  index  irregularities 
introduce  a  variety  of  signal  effects  including  multipath  and  group  delay 
jitter,  angle  of  arrival  jitter  and  spatial  decorrelation,  we  will  focus  on 
what  are  probably  the  most  important  effects,  namely  signal  amplitude  and 
phase  fluctuations. 

In  Fig.  4  we  illustrate  the  propagation  of  a  radio  signal  through 
an  irregular  medium  consisting  of  multiple  layers,  each  containing 
Irregularities  in  the  plasma  electron  number  density  N(r)  and  hence  in  the 
radio  refractive  index  n(r)  .  As  the  radio  ray  path  moves  relative  to  the 
Irregularities  the  received  radio  signal  scintillates,  i.e.,  undergoes 
variations  in  amplitude  and  phase  with  time.  The  twinkling  of  stars  is 
another  example  of  scintillation  caused  by  propagaticn  through  an  Irregular 
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Figure  4  Geometry  for  a  sir>fptified  propagation  model.  As  the  radio  ray  path  from  satellite 
to  Eardi  moves  relative  to  die  refractive  index  irregularities  (modelled  as  layers  in 
the  ionosphere),  these  irregularities  introduce  random  fluctuations  in  the  receivad 
signal's  amplitude  and  ph«e  (see  Figure  5)  known  as  signal  scintillation. 


A  useful  way  to  think  of  the  scintillation  process,  described  by 
Uscinskl  (1977,  Ch.  A),  Is  shown  In  Fig.  5.  Here  we  visualize  the  complex 
electric  field  £(2)  observed  at  some  point  2  along  the  straight  line  ray 
path  from  transmitter  to  receiver  (see  Fig.  4)  as  being  composed  of  a  non¬ 
scintillating  component  <£>  and  a  scintillating,  i.e.,  randomly 
fluctuating,  component  £  such  that 


£(2)  =  <fi(2)>  +  £  (2) 


(II-l) 


The  average  power  in  the  wave  field  is  therefore 


<££*>  «  <£>^  +  <££*>  . 


(I 1-2) 


At  the  top  of  the  scintillating  medium,  but  still  in  free  space,  we  take  z 

to  be  zero  increasing  downward  through  the  medium.  At  z  *  0  the  power  in 

2  2 

the  wave  field  in  all  in  the  coherent  component  <E>  ,  i.e.,  <EE  >  «  <E>  *  E 

As  z  increases  <£(z)>^  decays  exponentially: 


<£>^  »  E^  exp(-Bz) 


(11-3) 


where  can  be  thought  of  as  ’he  mean  free  path  for  radio  photons 

propagating  through  a  layer  of  the  scintillating  (irregular)  medium 

(Uscinskl,  1977,  ch.  4).  After  propagating  some  distance  ~  through 

the  mediinti  <£>^  has  decreased  by  e*"^  and  some  power  is  now  resident  in 

the  fluctuating  component  <£  £  >  as  shown  in  Fig.  5a.  Since  the 

variance  0?  of  the  phase  of  £  is  still  small  (o^  <  1  radian),  the 
9  9 
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Figure  S  Argand  (phasor  or^complex  plane)  diagram  ihowtng  the  tctai  electric  field  E  beiiw  composed 
of  a  meanfieM  <E>  along  the  abscissa  (real  axis)  and  a  fluctuating  component  E  such  that 
<E>  -fE'.  For  weak  scatur  (a)  amplitude  fluctuations  are  small  dll' l«i<^>  Hand 
phase  fluctuation  much  less  than  a  radian.  In  the  case  of  intermediate  scatter  (b)  amplitude 
fluctuations  arc  still  small,  but  phase  fluctuations  have  become  large.  For  strong  scatter  (c)  both 
amplitude  and  phase  fluctuations  are  large  (iTki  >  I);  signal  sutistics  approach  a  Rayleigh 
distribution  for  amplitude  and  uniform  for  phase. 


signal  Is  said  to  be  only  veakly  scattered.  Although  we  take  6  to  be 
constant  here  for  simplicity.  It  could  of  course  be  a  functi^'n  of  z  . 


A  variety  of  methods  exists  for  dealing  with  wave  propagation 
through  a  random  medium,  e«g.,  soms  recent  work  on  phase  screens  by  Rlno, 
Vondrak  and  Hatfield  (1975),  Rlno  (1979a  &  b)  and  Vlttwer  (1980a)  and  on 
Feynman  path  Integrals  by  Flatt4  et  al*  (1979) •  Of  primary  Importance  In 
all  these  methods  Is  the  spatial  size  distribution  of  the  Irregularities  in 
refractive  Index,  l.e.,  the  blobs  shown  In  Fig.  4.  For  weak  scatter 
fluctuations  In  amplitude  |  and  hence  In  power  <?£*>,  only 
refractive  index  fluctuations  on  spatial  scales  (£)  of  the  order  of  the 
Fresnel  zone  size  (In  the  scintillating  medium)  and  smaller  are 
Important.  In  the  case  of  Ionospheric  amplitude  scintillation  for  DHF 
signals  (X  ~  0.5  m)  only  size  scales  of  5  X  3  X  10^  «  400  m 

are  Important.  R  Is  the  distance  from  the  receiver  to  the  scattering 
medium.  As  scattering  becomes  strong,  larger  size  scales  become 
significant  for  amplitude  scintillation.  In  the  case  of  phase 
scintillation,  refractive  Index  fluctuations  on  size  scales  larger  as  well 
as  smaller  than  the  Fresnel  zone  size  are  Important  for  both  strong  and 
weak  scattering.  The  upper  limit  on  Important  size  scales  for  phase 
scintillation  Is  set  by  the  length  of  time  during  which  phase  is  observed » 
relative  velocity  of  the  ray  path  through  the  scintillating  mediiim  and  the 
propagation  geometry  (Rlno,  1979a)*  Since  fluctuations  of  significant 
strength  exist  down  to  size  scales  of  <  10  m,  a  range  of  size  scales  well 
over  an  order  of  magnitude  can  be  lii^>ortant  In  determining  the 
scintillation  of  radio  signals  propagating  through  a  disturbed  ionosphere. 

As  discussed  above,  a  signal  with  Raylei^  distributed  amplitude 
and  uniformly  distributed  pnase  (l.e.,  saturated  scintillation)  can  be  used 
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as  a  bounding  apprc ^Imatlon  for  signals  which  are  strongly  scattered  upon 
passage  through  a  disturbed  ionosphere  (see  Figs.  4  and  5).  Fig.  6  shows 
an  example  of  saturated  phase  and  amplitude  scintillation  on  a  UHF  signal 
propagating  through  a  disturbed  Ionosphere.  Communication  using  such  a 
Rayleigh  fading  channel  is  not  unique  to  transionospheric  satellite  links, 
but  is  encountered  in  communications  via  HF  skywave  as  well  as  VHF 
ionospheric  and  tropospheric  scatter  circuits.  Because  of  these  many 
practical  applications,  Rayleigh-fadlng  channels  have  been  studied 
extensively  (Schwartz,  Bsnnet  and  Stein,  1966).  Such  a  Raylelgh-fading 
channel  can  be  described  to  first  order  by  two  parameters  and  f^  . 

The  signal  decorrelation  time  measures  the  autocorrelation  property  of 

the  complex  electric  field  signal  at  the  receiving  antenna,  i.e., 
corresponds  to  the  time  lag  requin^d  for  this  autocorrelation  function  to 
drop  by  a  factor  e  from  Its  value  at  zero  time  lag.  Thus  x^  gives  the 
time  scale  over  which  the  E  vector  in  Fig.  5  moves  along  its  "ball  of 
string"  paths  in  the  Argand  diagrams.  In  Fig.  6  x^  corresponds  to  the 
mean  duration  of  fades  at  the  -3  dB  level.  The  other  parameter  f^  , 
called  the  frequency  selective  bandwidth,  measures  the  correlation  of 
fading  characteristics  at  two  separated  frequencies.  DNA  defines  f^  as 
l/2ito^  where  is  the  variance  of  the  random,  time-delay- jitter 

(excluding  jitter  associated  with  the  total  electron  content  along  the  ray 
path).  By  this  definition  f^  is  the  frequency  separation  required  to 
reduce  the  complex  signal  correlation  by  a  factor  0.697  below  its  value  at 
zero  frequency  separation.  These  parameters  are  the  primary  output  of 
computational  simulations  of  radio  signal  propagation  through  a  nuclear 
disturbed  ionosphere.  The  spatial  distribution  and  strength  of  plasma 


Figure  6  Example  of  UHF  signal  amplitude  and  phase  scintillation  for  strong  scatter  conditions  in  a 
dinurbed  ionosphere  with  signal  decontlation  time  «  8  ms.  Signal  statistics  approach 
the  saturated  regime  in  which  amplitude  is  Rayleigh  distributed  and  phase  is  uniformly 
distributed.  Scintillations  can  be  faster  or  much  slower,  depending  on  link  geometry  and 
ray  path  velocity  with  respect  to  the  plasma  irregularities.  After  McDaniel  (1978). 
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Irregularities,  ray  path  speed  with  respect  to  the  irregularities, 

propagation  path  geometry  and  operating  frequency  all  influence 

T  and  f 
o  o 

A  communications  channel  can  be  described  by 

R(t)  ■=  j  h(t,T)  S  (t  -  t)  dr  CII-S) 

o 

where  R  is  the  received  signa? ,  S  the  transmitted  signal  and  h  the 
impulse  response  of  the  channel.  Aside  from  signal  effects  associated  with 
propagation  geometry,  absorption.  Integrated  electron  content  along  the  ray 
path  and  angle  of  arrival  variations,  h(t,T)  for  the  Rayleigh-f ading  case 
can  be  characterized  by  the  two  parameters  and  f^  .  Since  for  most 

communications  systems  scintillation  effects  parametrized  by  and  f^ 

are  most  important  and  since  knowledge  of  h  enables  one  to  calculate 
error  ’■ates  for  various  communication  systems,  operating  frequencies, 
modulation  schemes,  etc.,  scintillation  effects  on  communications  systems 
can  be  at  least  bounded  by  knowing  x^^  and  f^  .  Further  discussion  is 
given  by  Wittwer  (1980a  &  b). 

In  Figs.  7  and  8  we  show  as  a  function  of  operating  frequency  the 
ranges  of  x^  and  f^  expected  for  communications  through  a  nuclear 
disturbed  ionosphere.  For  reference,  FLTSATCOM  operate- i  at  UHF;  Global 
Positioning  System  (GPS)  at  L  band;  Defense  Support  Program  (DSP)  at  S 
band;  WESTAR,  RADUGA  and  INTELSAT  TV  (and  V)  at  C  band  (4-8  GHz);  Defense 
Satellite  Communication  System,  (DSCS)  and  INTELSAT  V  at  X  band 
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EXPECTED  RANGE  OF  r  FOR 
SATURATED  (Rayleigh) 
SCINTILLATION 
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\  ESTIMATED  RANGE  OF  ^ 
V  CORRELATION  BANDWIDTHs 
\  IN  SATURATED  (Rayleigh!  v 
^SCINTILLATION  ' 


1  L  S 


FREQUENCY  —  GHz 


Figure  8  Estimated  range  of  signal  correlation  (or  frequency  selective)  bandwidths  for  strong 
scintillation  conditions  encountered  in  a  nucleat  environment.  The  minimum  values 
represent  very  severe  environments  and  are  less  likely  to  be  encountered  on  a  given  link. 
After  McDaniel  (1978). 
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(8-12.5  GHz)  and  the  planned  Strategic  Satellite  System  (S^)  at  EHF.  The 
frequencies  corresponding  to  these  various  bands  are  shown  In  Figs.  7  and 

8. 


Although  various  system  parameters  and  modulation  schemes 

determine  the  detailed  structure  of  error  rates  In  terms  of  t  and  f  , 

o  o 

some  qualitative  feeling  for  the  severity  of  scintillation  effects  can  be 
obtained  by  comparing  and  f^  with  the  time  scale  T  for  significant 
signal  variation.  T  would  usually  be  equivalent  to  the  time  required  to 
transmit  one  bit  of  Information  with  (1/T)  being  the  Information  trans¬ 
mission  rate  or  one  half  the  required  system  bandwidth  (except  for  spread 
spectrum  systems.)  Thus  for  all  systems  one  would  expect  catastrophic 
disruption  for  x  <  T  and/or  f  <  (2/T)  »  l.e.,  for  fades  lasting  on  the 

O  ^  O 

order  of  the  symbol  period  or  less  and  for  frequency  selective  bandwidths 

less  than  the  required  system  bandwidth.  For  coherent  systems  catastrophic 

degradation  occurs  for  <  1/Bj^  ,  where  Is  the  phase-tracking  loop 

bandwidth,  because  the  system  cannot  handle  the  rapid  phase  changes.  For 

T  >  T  (or  T  >  1/B,  for  coherent  systems)  amplitude  scintillation  will 
o  o  L 

still  cause  some  degradation.  For  example,  a  communication  system 
operating  at  one  megabit  data  rate  at  S  band  could  certainly  anticipate 
problems  In  view  of  Figs.  7  and  8  If  remedial  action  wer--.  not  taken.  By 
remedial  action  we  refer  to  various  mitigation  techn/'.p.is  discussed,  for 
example,  by  M'^Daniel  (1978)  and  Wlttwer  (1980a).  These  includi!  frequency 
shift  or  differential  phase  8h,(^  nvy'.~g  In  place  of  absolute  phase  shift 
keying,  high  transmitter  p;'.'. -.largln),  large  loop  bandw:.dth  phase 
trackers,  higher  operatl^ig  fx&quencles,  etc- 
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To  emphasize  the  importance  of  scintillation  effects  we  illustrate 
in  Fig.  9  the  large  area  of  the  United  States  which  would  be  affected  by  a 
high-altitude  nuclear  explosion  similar  to  that  shown  in  Fig.  1. 


Figure  9  Geographical  extent  of  strong  scintillation  along  links  between  ground  and  synchronous 
satellite  terminals  one  half  hour  after  a  high  altitude  nuclear  burst.  Note  the  wide  geo¬ 
graphic  area  affected  by  a  single  burst  and  the  strorrg  frequency  depertdence.  The  extent 
of  EHF  strong  scintillation  is  particularly  sensitive  to  uncertainties  arxl  nnay  be  smaller 
than  shown.  After  McDaniel  (1978). 


number  of  structures  such  as  shown  In  Fig.  2  and  after  scattering  from 
about  six  such  structures  the  resultant  wave  is  Indistinguishable  from  a 
wave  propagating  through  an  analagous  environment  described 
statistically.  Second,  because  one  would  not  know  just  exactly  \d\at  sort 
of  deterministic  structures  would  exist  in  a  nuclear  disturbed  ionosphere 
except  in  a  statistical  sense,  a  statistical  model  represents  "average" 
conditions. 


A  convenient  way  to  describe  the  distribution  of  the  strength  of 

refractive  index  irregularities  with  respect  to  size  is  to  use  a  spectral 

density  function  $  where  the  size  (6x)  of  an  irregularity  is  given  in 

terms  of  a  spatial  wavenumber  k  •=  2ti/6x  .  As  discussed  above,  refractive 

index  fluctuations  6n  in  the  ionosphere  arise  because  of  electron  density 

fluctuations  6N  where  for  small  values  of  (f  /f)  ,  6n  »  . 

From  this  we  see  that  the  ®  for  refractive  index  fluctuations  (6n)  is 

simply  related  to  the  ♦  for  electron  density  (or  plasma  density  in  the 

ionosphere)  fluctuations  6N  .  Further,  this  relation  is  scaled  by  the 

Inverse  fourth  power  of  the  operating  frequency  f  ,  i.e., 

$  «  (40.3)^  f~^  (inks).  Ionospheric  plasma  irregularities  arising  from 

n 

both  natural  causes  and  high  altitude  nuclear  explosions  are  reasonably 

well  fitted  by  power  law  k  °  curves  when  described  in  one-dimensional 

spatial  wavenumber  space,  i.e.,  by  measurements  along  some  straight  line  in 

space.  Thus  the  irregularities  can  be  described  by  #  (k;t,r)  where  ♦ 

is  the  spectral  density  function  for  refractive  Index  fluctuations,  k  Is 

♦ 

the  three  dimensional  wave  number  vector,  t  is  time,  r  gives  location 
along  the  ray  path  and  the  three-dimensional  spectral  index  (p)  is 
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related  to  the  one-dlmenslonal  spectral  Index  (a)  by  p  ■>  a  ■<-  2  for 

Isotropic  fluctuations  (Fejer  and  Kelley,  1980;  Rlno,  1978).  Since  the 

power  law  region  is  usually  bounded  by  a  lower  limit  ^  ibelow  which  the 

o 

wavenumber  spectrum  is  relatively  constant  %rlth  respect  to  jc  and  an  upper 
limit  above  which  the  spectrum  decreases  rapidly  with  Increasing  ^  , 

It  Is  reasonable  to  describe  such  plasma  irregularities  by  analogy  with 
fluid  turbulence  although  the  physical  processes  Involved  are  different. 

For  three^-dlmenslonal.  Isotropic  turbulence  we  can  generalize  the 
von  Karman  spectrum  (Ishlmaru,  1978,  ch.  16)  to  arbitrary  spectral  Index 
and  write 


♦(k;t,r) 


C^(k^  +k^)"^^“'^  ^^^^^exp(-k^/kj) 
o  i 


(II-6) 


where  C  is  a  "structore  constant”  describing  the  strength  of  the 
turbulence ;  k^  «  l/L^  anr.  is  the  outer  scale  size;  k^  •  1/^  and 

Is  the  Inner  scale  size;  and  a  Is  the  one~dlmensional  spectral 
Index.  In  terms  of  tvj'‘bulent  eddies  is  the  size  of  the  largest  eddy 

and  the  size  of  the  smallest.  In  Kolomogorov  fluid  turbulence 

a  ■>  5/3  .  For  plasma  striations  such  as  those  of  Fig.<  2  which  are 
dominated  by  sharp  gradients,  a  ~  2  .  A  sketch  of  (II-6)  is  given  In 


10  Sketch  of  the  turbulence  model  given  by  eg.  (11-6).  The  outer  scale 

corresponds  to  l/L^  while  the  inner  scale  Lj  corresponds  to  kj  “  1/Lj. 

For  Kolmogorov  fluid  turbulence  She  regions  labeled  1, 2  and  3  are  referred  to 
as  the  input  range,  the  inertial  subrange  and  the  dissipation  range  respectively 
and  a  -  5/3.  For  the  plasma  irregularities  considered  here  different  physical 
processes  are  at  work  and  0*^2. 


Due  to  the  geomagnetic  field  in  the  ionosphere,  irregularities  are 
not  isotropic  as  one  can  easily  see  in  Fig.  2.  Hence  one  must  generalize 
Eq.  (II-6)  to  include  scale  lengths  along  three  directions  as  done  by 
Wlttwer  (1980a); 


♦(k) 


3/2  2 

8;r  <An  >L  L  L  [r(0.5a  +  l)r(0.5o  -  0.5)] 
X  y  z 


,  2  2  2  2  2  2. 
(l+Lk  +Lk  +Lk) 
XX  y  y  2  z 


(o+  2)/2 


(II-7) 


2 

where  <An  >  is  the  total  variance  of  the  refractive  index  fluctuations; 

,  Ly  ,  and  are  the  outer  scales  sizes  in  the  x  ,  y  ,  and  z 

directions,  F  denotes  the  Gamma  function  and  a  the  one-dimensional 
spectral  index.  Although  no  inner  scale  is  included  in  Eq.  (11-7)  one 
could  be  included  as  in  Eq.  II-6  by  using  a  multiplicative  exponential 
term. 

Once  the  spatial  wavenumber  spectrum  of  refractive  index 
irregularities  is  specified  as  in  equations  11-6  and  11-7,  the  spatial  and 
temporal  characteristics  of  signal  amplitude  (or  intensity)  and  phase 
fluctuations  can  be  calculated.  That  is,  one  can  calculate  the  signal 
decorrelation  time  and  frequency  selective  bandwidth  f^  as  well  as 

many  other  signal  statistics,  for  example,  the  temporal  frequency  spectra 
of  signal  amplitude  and  phase  fluctuations.  For  phase  fluctuation 
statistics  the  detrend  time  or  length  of  observations  in  time,  whichever  is 
shorter,  is  also  Important  (Rino,  1979a). 
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It  is  instructive  to  consider  the  simple  case  of  Fig.  11  where  a 

plane  wave  propagates  through  a  statistically  homogeneous  and  Isotropic 

medium.  The  Irregular  medium  and  ray  path  are  moving  relative  to  each 

other  with  a  velocity  which  is  transverse  to  the  ray  path.  In  this 

simple  case  of  weak  scattering  the  temporal  frequency  spectra  of  amplitude 

P.((jj)  or  phase  P,(a>)  fluctuations  can  be  written  as  integrals  of  the 
A  (p 

wavenumber  spectrum  of  refractive  index  fluctuations  tCk)  weighted  by  the 
appropriate  filter  function.  From  Ishimaru  (1977,  ch.  19) 


Pa  = 
A,<t) 


8Tt^<n>^a)^L 


c^V. 


f~f.  ,(k)*(k)dk'  (II-8) 

J  A,  4. 


where  <n>  is  the  average  refractive  index,  u  the  wave  frequency,  L 
the  length  of  the  ray  path  (the  scintillating  medium  covers  the  length  L 
of  the  ray  path),  c  the  velocity  of  light,  the  velocity  of  the  ray 

path  with  respect  to  the  medium  and  k  =  |^k'^  +  •  'P*'®  filter 

functions  f  and  f  for  amplitude  and  phase  are  illustrated  in  Fig.  12. 

A  (p 


The  important  feature  of  the  filter  functions  is  that  for  weak 
amplitude  scintillation  the  wavenumber  spectrum  of  irregularities  is 
effectively  high-pass  filtered,  l.e.,  the  influence  of  ♦  at 
k  <  ('\/2ir/ )  is  restricted.  Put  another  way,  only  refractive  Index 
fluctuations  at  scale  sizes  less  than  the  Fresnel  zone  size  are 


effective  in  producing  signal  amplitude  fluctuations  for  the  weak 
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Figure  11  Simple  example  of  plane  wave  propagating  through  a  statistically 
homogeneous  and  isotropic  scintillating  madium.  The  ray  path 
(and  hence  transminer  and  receiver)  move  with  a  transverse 
velocity  Vf  relative  to  the  medium. 


Vu  Vu 


Figure  12  General  shapes  of  weak  scatter  weighting  or  filter  functions  f^  and  f^  for 

signal  amplitude  and  phase  respectively.  Note  that  in  performing  the  integral 

of  eg.  II-8  the  functions  f  and  4>  are  evaluated  at  k  such  that  k^  >  k*^  -f 
2  2 

(u;  /V^^)  where  k'  is  the  variable  of  integration.  Figure  taken  from  Ishimaru 
(1977.  ch.  17). 
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scattering  case*  For  strong  scatter^  size  scales  larger  than  Fresnel  zone 

size  become  involved  in  amplitude  scintillation.  On  the  other  hand  all 

scale  sizes  potentially  Influence  signal  phase  fluctuations  in  either  weak 

or  strong  scatter.  However,  Rino  (1979a)  has  pointed  out  that  the  largest 

scale  sizes  which  can  influence  phase  fluctuations  are  limited  by  the 

detrend  time  or  total  observation  time,  whichever  is  shorter.  The  point  we 

make  here  is  that  the  wave-nu  ■•.ber  spectnnn  of  refractive  index  fluctuations 

4(k)  along  with  various  known  geometrical  and  radio  system  parameters  are 

the  essential  Ingredients  used  in  calculating  the  signal  statistics  one 

would  expect  upon  propagation  through  an  ionosphere  disturbed  by  a  high 

altitude  nuclear  explosion  (or  by  natural  ionospheric  disturbances,  for 

that  matter).  Although  we  have  illustrated  this  point  by  discussing 

P.  ^(ui)  ,  it  applies  equally  well  to  other  signal  statistics  such  as 
A,(j> 

and  f^  ,  the  signal  decorrelation  time  and  frequency  selective 
bandwidth.  Since  the  propagation  geometry  and  radio  system  factors  are 
known  and  slice  the  radio  propagation  theory  is  relatively  well  understood, 
the  weak  link  in  the  calculation  turns  out  to  be  ♦(k)  ,  the  spatial 
distribution  of  refractive  index  irregularities.  Hence  the  bulk  of  this 
report  focusses  on  the  determination  of  4(k)  ,  the  ultimate  objective 
being  to  accurately  model  4(k)  as  a  function  of  time,  spatial  location 
and  scenario. 

C.  Sensitivity  of  Signal  Statistics  to  Plasma  Irregularity 
Structure 

Now  that  we  have  identified  the  spatial  wavenumber  spectrum 
of  refractive  index  irregularities  4(k)  as  a  weak  link  in  our  ability  to 


11-23 


calculate  radio  signal  propagation  through  a  disturbed  ionosphere,  it  is 
appropriate  to  ask  how  sensitive  signal  statistics  are  to  errors  in  our 
knowledge  of  4'(k)  ,  i.e.,  to  errors  In  the  turbulence  strength,  outer 
scale  Lq  and  spectral  index  a  •  Wlttwer  (1980c)  has  calculated  signal 
decorrelation  time  (t^)  and  frequency  selective  (l«e.,  available) 
bandwidth  (f^)  as  functions  of  the  spectral  index  a  and  outer  scale 
size  Lq  •  These  calculations  are  for  a  nominal  Rayleigh  fading  case. 

When  is  the  independent  variable,  a  is  held  at  2.  Both  and 

Iq  can  be  closely  approximated  by  power  law  functions  of  over  the 

range  of  encountered  in  the  ionosphere,  namely 

“  [L^/^na^/L^))®*^  or  i=-  1.7  a  L^/toCL^/L^)  or 

f  =1.3  for  *  10  m  and  1  km  <  L  <  30  km  where  x  is  in 

msec,  fg  in  MHz  and  Lq  in  km.  Both  x^  and  f^  turn  out  to  be  strong 
functions  of  the  spectral  index  a  as  illustrated  in  Fig.  131 


It  is  clear  from  these  nominal  case  calculations  that  signal 
statistics  are  indeed  sensitive  to  errors  in  the  parameters  characterizing 
the  spatial  distribution  of  refractive  index  irregularities  in  the 
disturbed  ionosphere  4>(?)  .  The  spectral  index  a  has  a  particularly 
strong  influence  for  a  <  2  .  The  current  canonical  value  for  o  is  2 
for  late  times.  If  plasma  irregularities  were  to  evolve  at  late  times  such 
that  a  is  substantially  below  2  ,  then  current  projections  which  take 
a  »  2  would  be  seriously  underestimating  the  signal  distortion  produced 
by  propagation  through  a  disturbed  ionosphere.  The  evolution  of  and 
a  with  time  are  discussed  in  Section  III  below. 
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SIGNAL  DECORRELATION  TIME  (r„) 


SPECTRAL  INDEX  (o) 


Figure  13  Signal  decorrelation  time  (Tg,  solid  curve)  and  frequency  selective  bandwidth  (f^, 
dashed  curve)  as  functions  of  the  spectral  index  Q  of  a  power  law  model  for  the 
wavenumber  spectrum  $(lc)  of  refractive  index  fluctuations.  After  Wittwer  (1960c). 
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D.  Effects  of  Plasma  Irregularities  on  HF  Communications  via 
Ionospheric  Refraction  and  Scattering 

1 .  Importance  of  HF  Communications  Effects 

Although  satellite,  cable  and  microwave  relay  links 
carry  the  bulk  of  military  and  civil  communications,  point  to  point 
communications  by  HF  (~  3  to  30  MHz)  sky  wave  (ionospheric  refraction) 
propagation  are  frequently  the  best  (and  often  the  only)  back-up  system 
available  for  long  distances.  Soviet  forces  apparently  rely  on  HF  systems 
for  much  of  their  primary,  as  well  as  secondary,  communications  needs.  If 
HF  sky  wave  systems  are  to  provide  viable  back-up  communications  for  use 
when  satellite  links  are  unavailable,  they  must  naturally  function  as  well 
as  possible  when  the  ionosphere  is  disturbed  by  high  altitude  nuclear 
explosions.  Many  nuclear  effects  on  HF  communications  have  been 
extensively  studied,  e.g.,  increased  absorption  due  to  increased  ionization 
at  relatively  low  altitudes.  Propagation  codes  which  model  such  large 
scale  effects  have  existed  for  over  a  decade,  e.g.,  NUCOM  (Owen  and  Bailey, 
1978;  Baumann,  Hatfield  and  Owen,  1973;  and  Neilson,  Lomax  and  Turner, 
1967).  However,  to  the  author's  knowledge  such  codes  do  not  include  the 
effects  of  scattering  from  small  scale  (<  50  km)  irregularities.  Since 
such  small  scale  ionospheric  irregularities  will  presumably  exist  for  hours 
after  a  high  altitude  nuclear  explosion,  it  is  necessary  to  assess  their 
impact  on  existing  and  planned  HF  sky  wave  communications  systems  which 
must  operate  in  a  nuclear  environment.  Below  we  discuss  the  probable 
nature  of  these  scattering  effects,  an  approach  to  dealing  with  them 
theoretically  and  how  they  might  be  used  in  a  VHF  (35-50  MHz)  ionospheric- 


scatter  communications  system. 


2.  The  Impact  of  Scattering  by  Small -scale  Ionospheric 
Irregularities  on  HF  Sky-Wave  Communications 

Since  Ionospheric  irregularities  caused  by  high  altitude 

nuclear  explosions  are  in  many  ways  analagous  to  ionospheric  irregularities 

occurring  naturally  in  equatorial  and  auroral  regions ,  one  can  argue  by 

analogy  that  study  of  sky  wave  communications  circuits  propagating  through 

these  natural  irregularities  should  provide  information  relevant  to  the 

nuclear  case*  The  ionospheric  irregularities  are  identified  by  a  condition 

known  as  "spread  F."  Spread  F  refers  to  ionospheric  conditions  in  which 

vertically  transmitted  HF  waves  are  reflected  not  by  a  single  well  defined 

ionospheric  layer,  but  by  irregularities  spread  over  a  range  i’.  altitude* 


Further  information  on  spread  F  and  ionospheric  irregularities  is  given  by 
Davies  (1965)  and  Fejer  and  Kelley  (1980)* 


HF  sky  wave  propagation  paths  passing  through  spread  F  regions, 
l.e*,  regions  containing  ionospheric  irregularities,  are  not  refracted  by 
the  ionosphere  in  a  smooth  fashion,  but  rather  a  number  of  wavelets 
originating  in  various  regions  arrive  at  the  receiving  antenna  where  their 
vector  sum  is  sensed  by  the  receiving  system*  Interference  between  the 
arriving  wavelets  causes  signal  amplitude  and  phase  to  fluctuate  in  a 
random  manner  as  with  the  satellite  signals  discussed  above*  In  the  sky 
wave  case  the  propagation  ray  path  is  relatively  stationarj'  and  relative 
movement  between  ray  path  and  ionospheric  Irregularities  is  caused 
primarily  by  Irregularity  movement* 


Experience  with  sky  wave  cofimunlcations  along  paths  where  spread  F 
is  common,  l.e*,  through  the  auroral  zone  or  through  equatorial  regions 
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(particularly  transequatorlal  paths),  shows  that  flutter  fading  Is  comon 
with  fading  rates  from  a  few  to  perhaps  100  Hz  (Davies,  1965,  Ch.  5). 

These  fading  rates  correspond  to  Irregularity  motion  of  ~  100  ms“^.  The 
frequency  selective  bandwidth  typically  ranges  from  ~  3  kHz  for  CW  signals 
to  ~  40  kHz  for  short  pulses  (op.  clt.).  Such  flutter  fading  Is 
disastrous  to  conventional  voice  communication. 

A  potential  source  of  further  experimental  Information  would  be  HF 
sky  wave  transmission  experiments  carried  out  during  the  1962  Fishbowl 
series  of  high  altitude  nuclear  tests.  Oblique  Ionospheric  sounders 
between  various  points  in  the  Pacific  were  run  for  periods  (~  hours)  after 
high  altitude  explosions  and  spread  F  type  conditions  were  observed,  e.g., 
along  the  2500  km  path  between  Canton  Is.  (2°  S,  172°  W)  and  Raratonga  Is. 
(21°  S,  159°  W)  from  about  50  minutes  after  the  Kingfish  test  to  130 
minutes  and  presumably  longer  (Nielson,  1972). 

One  should  bear  in  mind  here  that  fading  is  a  relatively  common 
phenomenon  on  HF  sky  wave  paths  and  arises  from  a  number  of  mechanisms  in 
addition  to  the  small-scale  irregularities  at  Issue  here.  Hence  the  users 
of  HF  sky  wave  communications  circuits  have  always  had  to  cope  with  signal 
fading  and  have  adjusted  system  designs  accordingly.  On  the  other  hand  the 
flutter  fading  on  auroral  and  transequatorlal  paths  is  very  severe  relative 
to  that  typically  experienced  at  mid-latitudes.  In  addition  fading  effects 
associated  with  a  nuclear  disturbed  ionosphere  are  likely  to  be  worse 
still.  In  sum  it  would  be  prudent  to  investigate  the  impact  of  small-scale 
ionospheric  irregularities  on  HF  sky  wave  circuits  which  are  to  operate  in 
a  nuclear  environment. 
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3.  A  Theoretical  Approach  to  Scattering  of  HF  Signals  by 

Small-Scale  Ionospheric  Irregularities 

The  effects  of  small-scale  irregularities  on  HF  wave 
propagation  in  the  ionosphere  were  studied  by  Pitteway  (1960a  &  b)  and  the 
results  applied  to  measurements  of  the  movements  of  ionospheric 
irregularities,  Wright  (1968  and  following  papers)  and  Pitteway  (1971). 
While  extension  and  application  of  Pitteway's  work  could  be  fruitful,  it  is 
probably  better  to  pursue  a  different  course,  taking  advantage  of  more 
recent  work  by  Flatt4  et  al.  (1979)  and  Haugstad  (1978).  The  problem  of  HF 
sky  wave  propagation  through  an  Irregular  medium  has  much  in  common  with 
the  translonospherlc  propagation  problem  discussed  above.  The  principal 
difference  is  that  in  the  HF  case  the  background  medium  (in  which  the 
irregularities  exist)  is  Important — it  refracts  the  HF  wave  back  to 
Earth.  In  the  translonospherlc  (satellite)  case  the  background  medium  is 
not  Important  in  this  sense  since  little  large-scale  refraction  occurs. 

The  Importance  of  the  background  medium  is  clear  in  the  case  of  sound 
transmission  in  Che  ocean  (SOFAR)  channel  and  hence  is  included  in  the  work 
of  Flatt4  et  al.  (1979).  Similarly  Haugstad  also  included  background 
effects  in  his  study  of  EM  wave  propagation  through  planetary 
atmospheres.  It  would  appear  that  HF  sky  wave  propagation  through  a 
disturbed  ionosphere  could  in  some  cases  be  treated  by  applying  the  work  of 
Flatt^  et  al.  and/or  Haugstad  to  a  model  of  the  ionosphere  disturbed  by  a 
nuclear  explosion  as  discussed  briefly  in  Section  I  above.  The  cases  where 
this  method  would  likely  fall  would  be  those  where  the  irregularities  %»ere 
very  strong,  causing  large  angle  scattering.  At  late  times  when  small 
angle  scattering  would  presumably  dominate,  these  methods  should  work. 
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A.  Ionospheric  Scatter  Communication  Using  Irregularities 

Produced  by  High  Altitude  Nuclear  Explosions 

Irregularities  in  the  Ionosphere  can  be  used  to 
communicate  at  30-50  MHz  over  ~  1500  km  paths  in  a  manner  similar  to 
tropospheric  scatter  systems  at  UHF  frequencies.  In  fact  an  ionospheric 
scatter  communications  system  was  constructed  by  the  Air  Force  stretching 
from  Massachusetts  via  Canada,  Greenland,  Iceland  and  Britain  to  Eurc^e 
(JTAC,  1960).  This  system  depended  on  relatively  large  antennas  (~  20  dB 
gain),  high  powers  (~  tens  of  kW)  and  allowed  communications  bandwidths 
relatively  small  (~  10  kHz)  compared  to  tropospheric  scatter  systems. 
Irregularities  at  D-region  altitudes  provided  the  scattering  medium.  The 
principal  advantages  over  troposcatter  lay  in  the  longer  path  lengths 
possible  and  more  simple  apparatus  required. 

The  ionospheric  irregularities  caused  by  high  altitude  nuclear 
explosions  would  provide  an  abundance  of  scatterers  which  could  be  used  for 
communications.  The  point  here  is  that  ionospheric  scatter  communications 
systems  are  practical  and  might  provide  a  viable  back-up  (with  respect  to 
satellites  and  landline)  mode  for  long  distance  (~  1500  km)  point  to  point 
communications  in  a  nuclear  environment.  While  the  data  rate  of 
ionospheric  scatter  channels  is  only  comparable  to  HF  sky  wave  channels, 
the  higher  operating  frequencies  would  reduce  absorption  along  the  path  due 
either  to  natural  phenomena  or  to  the  nuclear  environment.  In  adaltlon  the 
nuclear  explosion  related  Irregularities  at  F-region  altitudes  would 
presumably  provide  lower  loss  and  longer  distance  paths  than  the  E-region 
irregularities  used  in  the  natural  environment. 
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In  a  nuclear  disturbed  ionosphere  adaptive  HF-VHF  systems  might  be 
able  to  use  refractive  (as  opposed  to  scatter)  paths.  This  could  be  a 
better  option  than  the  scatter  system  suggested  above.  However,  neither 
option  has  been  explored  quantitatively  to  the  author's  knowledge. 
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III.  PHYSICS  OF  PLASMA  IRREGULARITY 
EVOLUTION  AT  LATE  TIMES 


A.  Important  Physical  Mechanisms  for  Late  Time  Evolution  of 
Striated  Plasmas 

We  are  concerned  here  with  plasmas  (whether  from  a  high- 
altitude  nuclear  explosion,  a  barium  release,  or  a  laboratory  simulation) 
whose  size  and  shape  scale  lengths  are  very  much  larger  than  a  set  of  scale 
lengths  identified  below.  Under  these  circumstances,  the  plasma  clouds 
evolve,  under  the  Influence  of  a  differential  electron-ion  drift,  by 
spontaneously  generating  smaller  scale  lengths.  This  process — usually 
idealized  as  the  ^  x  S  gradient-drift  instability  (Linson  and  Workman, 
1970;  Perkins  et  al.  1973,  and  Volk  and  Haerendel,  1971) — has  one  Important 
feature  in  common  with  inertial-range  turbulence:  the  evolution  equations 
make  no  reference  to  any  of  the  small  length  scales  of  the  problem.  We 
sketch  a  plausible  outcome  of  the  long-term  behavior  of  scale-free 
evolution  in  Sec.  Ill  B;  here  we  are  concerned  with  the  questions  of  what 
are  the  small  scale  lengths,  and  how  do  they  affect  the  evolution  from 
large  to  small  length  scales. 

In  general,  scale-free  evolution  yields  a  power-law  spectrum  for 
plasma  irregularities  [Eqs.  (11-6  &  7)  and  Sec.  Ill  B],  and  at  k-values 
comparable  to  the  Inverse  of  the  smallest  scale  length  the  spectrum  is  cut 
off  (see  Fig.  10).  It  is  not  known  whether  the  cutoff  is  exponential,  as 
in  Eq.  (II-6),  but  it  hardly  matters,  as  long  as  the  cutoff  is  sharp 
enough.  One  does  not  knew  a  priori  exactly  how  the  cutoff  in  k-space 


depends  on  the  inverse  scale  length,  and  large  ratios  can  occur  (like 
turbulent  Reynolds  numbers),  so  that  one  must  be  prepared  to  study  the 
dynamics  associated  with  a  given  physical  scale  length. 


In  Table  I,  we  list  candidate  scale  lengths  for  the  three  types  of 
plasma  clouds  of  Interest  to  us  at  times  well  after  cloud  formation.  These 
length  scales  are  found  by  simple  dynamical  arguments,  e.g. ,  finding  a 
damping  terra  which  just  balances  the  linear  x  S  gradient-drift  growth 
rate.  We  comment  briefly  on  the  effects  listed  in  the  Table. 


Electron  Diffusion.  To  our  knowledge,  this  is  the  only  physical 
scale  length  whose  effects  are  claimed  to  be  observed  in  large-scale 
simulations  (McDonald  et  al.  1978).  Roughly  speaking,  the  electron 
diffusion  coefficient  Dg  is  given  by 


Ve  (VTj) 

°e  ~  2  m 

Q  e 


(III-l) 


e 

where  v  is  a  collision  frequency,  fi  a  gyrofrequency ,  T  a  temperature 
(in  energy  units),  and  e,  i  refer  to  electrons  and  ions.  To  find  the 
scale  length  equate  k^D  to  the  linear  growth  rate  y  =  kj^U  (D  ®  cE^/b, 
the  convective  drift  velocity).  It  is  easy  to  superimpose  diffusion 
dynamics  on  the  scale-free  equations,  because  neither  kinetic  effects  nor 
three-dimensional  geometries  need  be  considered,  and  presumably  this  is  the 
reason  why  electron  diffusion  is  the  only  process  Included  in  the  large 
computer  simulations.  This  is  an  unfortunate  state  of  affairs,  because 
electron  diffusion  is  the  smallest  scale  length  listed  in  the  Table  and  one 
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expects  Important  phenomena  at  the  larger  scale  sizes  listed  there.  [Inner 
scale  lengths,  which  are  substantially  smaller  than  the  Fresnel  zone  size 
of  ~  400  m  are  not  likely  to  be  Important  so  far  as  radio  communication 
systems  are  concerned.] 

Ion  Diffusion.  The  analog  of  (IIl-l)  for  ions  yields  a  diffusion 
coefficient  (and  scale  length)  about  ten  times  as  large,  but  it  is 
inappropriate  to  invoke  this  length  scale  for  a  two-dimensional  problem. 

The  well-known  reason  is  that  ions  cannot  diffuse  any  faster  than  the 
electrons  can  move  across  field  lines;  that  is,  normally  ion  diffusion  is 
exactly  cancelled  by  the  effect  of  the  polarization  potential  in  the 
cloud.  However,  if  the  plasma  cloud  is  immersed  in  a  background  (e.g., 
ionospheric)  plasma  of  sufficiently  high  conductivity,  electrons  can  in 
effect  move  across  field  lines  by  short-circuiting  through  the  background 
plasma  (Volk  and  Haerendel,  1971,  and  Francis  and  Perkins,  1975).  This  has 
been  described  in  detail  by  Francis  and  Perkins  who  show  that  if  the 
height-integrated  background  Pedersen  conductivity  is  large  enough, 

ion  diffusion  is  allowed  but  at  a  rate  reduced  by  the  factor 

L,  (Z,  +  Z  )”^  (  Z  is  the  plasma  cloud  integrated  conductivity).  On  the 
b  b  c  c 

other  hand,  these  authors  show  that  a  relatively  compressible  (collisional) 

background  may  actually  enhance  strlatlon  growth,  rather  than  lead  to 

diffusive  damping.  For  nuclear  disturbances  of  the  ionosphere  Z^  is 

usually  much  j.  irger  than  Z.  with  £  becoming  relatively  smaller  at  late 

b  c 

times. 


III-A 


Perkins,  Zabusky,  and  Doles  (1973)  have  described  another  somewhat 
similar  mechanism  involving  parallel  resistivity:  the  electrostatic 
potential  is  not  constant  along  a  magnetic  field  line.  In  this  case,  one 
loses  the  constraint  that  all  electrons  on  a  given  field  line  have  the  same 
X  drift,  and  limited  transverse  ion  diffusion  can  occur.  The  ion 
diffusion  scale  length  for  barium  clouds  given  in  Table  I  is  based  on  the 
discussion  in  Perkins  et  al.  (1973).  We  do  not  give  an  ion-diffusion 
length  for  the  nuclear  or  laboratory  cloud  in  the  Table,  because  of  many 
uncertainties  in  the  properties  of  the  background  plasmas  which  determine 
whether  ion  diffusion  is  important  or  not. 

The  important  fact  about  ion  diffusion  is  that  it  is  governed  by 
three-dimensional  effects,  which  clearly  complicate  large-scale  computer 
simulations.  There  have  been  some  attempts  to  mimic  three-dimensional 
effects  by  doing  two-level  simulations  [a  uniform  background  cloud  below  a 
uniform  plasma  cloud,  (Scannapieco  et  al.  1976)],  but  we  are  not  aware  of 
any  in  which  ion  diffusive  damping  could  have  been  seen. 

Ion  Inertia  Effects.  Often  one  is  concerned  only  with  cases  where 
is  the  largest  frequency  in  the  problem;  in  particular,  » 

where  v^  is  a  typical  ion  velocity.  Linear  ion-inertia  effects  are 
easily  accounted  for  (Linson  and  Workman,  1970,  and  Ossakow  et  al.  '978), 
but  non-linear  effects  which  might  affect  the  evolution  of  situations  whose 
size  comparable  to  the  Larraor  radius  R  *  ^  have  not  been 

considered  (as  far  as  we  know)  for  the  striation  problem,  although  the 
general  physical  picture  has  been  appreciated  for  very  similar  spread-F 


phenomena  (see  the  review  of  Fejer  and  Kelley,  1980;  also  Costa  and  Kelley, 
1978a  &  b). 

A  particularly  important  effect  is  the  coupling  of  the  gradient 
t  X  t  drift  instability  to  drift  waves.  The  scale  frequency  for  drift 
waves  is  oj*  : 


uS*  = 


mod 

i  i 


(111-2) 


where  d  is  a  typical  cloud  gradient  scale  length,  and  ~  1  for 

fastest  growing  modes  (if  the  drift  wave  is  unstable).  The  drift-wave 
frequency  is  comparable  to  the  gradient-^  x  S  growth  rate 
f  a  l]d“^  (U  =  cE^B"*^)  when  the  average  ion  velocity  is  comparable  to  U  , 
a  not  untypical  condition  for  the  plasma  clouds  of  interest.  In  order  that 
the  usual  linear  electrostatic  drift-wave  theory  apply,  for  which 
w*  <  ,  one  should  have  <  d  ,  but  it  will  be  necessary  to  consider 

the  evolution  of  gradient  scales  all  the  way  down  to  the  Larraor  radius  in  a 
complete  theory. 


In  the  absence  of  exotic  ion  distributions  in  velocity  space, 

drift-wave  instability  is  controlled  through  a  three-dimensional  effect;  a 

finite  k,,  such  that  k  v^  «  to*  <  k. ,v  .  This  means  that  the  cloud 
'  '  Hi  He 

must  be  no  longer  (along  the  field)  than  about  dv  v  »  d  ,  a  condition 

e  1 

which  may  or  may  not  be  satisfied. 


III-6 


It  Is  not  particularly  i'l^ortant  whether  these  modes  are  stable  or 


unstable.  Unstable  modes  will  eventually  be  damped,  by  a  variety  of 
mechanisms:  Ion  diffusion  (remember,  kj j  ^0)  ,  or  Interesting  non-llnear 
Interactions  like  coupling  to  convective  cells  (Okuda  and  Dawson,  1973) 
which  Is  a  kind  of  enhanced  diffusion,  or  coupling  to  what  Dupree  (1978) 
calls  clumps  (of  particles  In  a  BGK-llke  mode).  All  of  these  effects  can 
only  be  treated  adequately  by  taking  Into  account  some  features  of  kinetic 
theory,  which  raises  new  problems  for  computer  simulations.  (These 
problems  are  sufficiently  Important  and  specialized  that  we  have  devoted 
another  subsection  to  them.)  But  It  may  very  well  be  possible  to  model 
most  of  the  Interesting  drift-wave  effects  In  a  way  which  Is  amenable  to 
computer  codes  now  In  use  (2-d  electrostatic).  We  have  In  mind  a 
development  analogous  to  the  quasx-three-dlmenslonal  treatment  of  Hasegawa 
and  Mima  (1978),  who  take  into  account  the  parallel  motion  of  electrons 
along  the  magnetic  field,  yet  derive  completely  two-dimensional  non-linear 
electros; atic  equations  which  describe  important  features  of  drift-wave  and 


convective-cell  turbulence. 


We  believe  that  the  phenomena  described  here  under  the  heading  of 
ion  inertia  may  well  be  dominant  in  setting  the  Inner  scale  length  where 
the  gradient-associated  free  energy  of  the  E  x  B  gradient  drift 
instability  Is  dissipated. 


Higher-Frequency  Unstable  Modes.  There  is  a  very  considerable 


literature  on  such  modes  as  applied  to  F-region  irregularities  and  barium 


clouds  (Fejer  and  Kelley,  1980;  Huba  et  al.  1978;  Huba  and  Ossakow,  1979, 
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and  Simons  et  al.  1980).  We  will  not  go  into  details;  the  modes  invoked 

are  (besides  the  drift  waves  discussed  above)  the  electron  or  ion- 

cyclotron-drift,  lower-hybrid  drift,  and  drift-loss-cone.  Frequencies  are 

comparable  to  £1^  ,  and  wave  lengths  lun  from  the  ion  Larmor  radius  to  the 

electron  Larmor  radius.  The  theory  is  fully  kinetic,  thus  not  amenable 

(without  drastic  simplification)  to  the  usual  electrostatic  fluid 

simulation  codes.  (It  should  be  noted  that  some  of  these  modes  are  fully 

two-dimensional,  and  unstable  when  k| ^ =  0  .)  Whatever  effects  Lnese  modes 

have  can  be  accounted  for  in  particle-kinetic  simulations,  although  even 

these  simulations  will  be  difficult  for  electron-driven  modes.  i 


The  Los  Alamos  group  (Simonds  et  al.  1980)  has  interpreted  the 
prompt  striations  observed  in  the  barium  shaped-charge  cloud  Buaro  as  due 
to  a  drift-loss-cone  instability.  These  authors  have  suggested  to  us  in 
briefings  that  very  early-time  spatial  structures  of  size  <  ,  spread 

through  the  initial  cloud,  are  important  in  determining  the  way  the  late¬ 
time  gradient-fi  x  8  drift  instability  evolves.  However,  kinetic 
instabilities  such  as  considered  by  the  Los  Alamos  group  saturate  at  a 
modest  level  (6n/n  a  few  percent),  and  it  is  hard  to  see  how  these  are 
involved  in  the  100%  spatial  modulations  actually  observed  in  barium 
clouds. 


B.  The  Problem  of  Large  Dynamic  Range  and  Fractalization 

We  have  already  said  in  Sec.  1  that  the  plasma  physics  of  a 
high-altitude  nuclear  cloud  (or,  to  a  lesser  extent,  that  of  a  barium 
cloud)  covers  an  enormous  range  of  spatial  scales  at  late  times  (~  hours 
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after  cloud  f ormati' ■') .  This  range  is  so  large  (perhaps  up  to  10^)  that 
conventional  approaches  to  con^Juter  sitnulati'^n,  or  ev'm  the  largest 
laboratory  plasma  devices  (see  Sec.  IV),  cannot  begin  to  cover  it.  Just 
how  large  the  dynamical  range  is  depends  on  what  one  believes  governs  the 
physics  at  short  distances  (at  the  inner  scale);  Fig.  14  shows 
schematically  the  latc-time  divergence  of  the  outer  and  inner  scales  of  a 


nuclear  plasma  cloud  for  two  candidate  inner  scales:  the  ion  Larmon  radius 


and  the  electron  diffusion  length  (gotten  by  equating  k 


(y(?f 


the  gradient-drift  growth  rate  kj^U  ).  Also  shown  is  an  estimate  of  a 
computer's  dynamical  range,  based  on  a  100  x  100  spatial  mesh  and  assuming 


that  3  mesh  points  aie  needed  to  model  the  smallest  spatial  features. 


There  are  other,  larger  candidates  for  an  inner  scale  (discussed 
above),  and  a  barium  cloud  might  only  reach  an  outer  scale  of  10  km, 
instead  of  1000  km  for  a  nuclear  cloud,  so  there  are  circumstances  when  the 
dynamical  range  of  spatial  scales  for  an  ionospheric  barium  cloud  may  be 
accessible  to  straightforward  computer  simulation  or  even  laboratory  analog 
experiments.  But  it  is  clearly  of  interest  to  study,  in  a  semi-analytical 
way,  the  problems  raised  by  a  virtually  infinite  dynamic  range.  In  fact. 

It  is  common  in  the  computer  simulations  to  make  the  range  actually 
infinite  by  throwing  out  of  the  equations  of  motion  all  terms  which  refer 
to  a  fundamental  length  scale.  Since  the  resulting  scale-free  equations — 
termed  by  us  the  canonical  equations — have  only  a  few  specialized  analytic 
solutions,  there  has  been  little  progress  in  analysis  of  the  full  non¬ 
linear  problem. 


The  last  sentence  is  (or  was,  until  recently)  true  of  many 
physical  problems  in  quantum  field  theory,  hydrodynamlcal  turbulence, 
second-order  phase  transitions,  etc.  The  conmon  feature  is  an  enormous 
dynamical  range;  the  common  method  of  attack  is  to  exploit  the  self¬ 
similarity  of  the  scale-free  equations.  (That  is,  the  contours  in  Figs.  15 
and  16  are  all  made  of  figures  of  the  same  shape,  but  sizes  which  get 
smaller  as  the  fractal  evolves.  A  fractal,  discussed  more  fully  below,  is 
a  continuous  but  non-dif ferentiable  curve  of  infinite  length.) 

We  have  carried  out  a  preliminary  investigation  of  the  self- 
similarity  properties  of  the  canonical  equations,  with  the  result  that  a 
plasma  cloud  subject  to  these  equations  fractalizes,  that  is,  isodensity 
contours  (always  enclosing  the  same  area)  bifurcate  to  infinitely  small 
spatial  scales  and  become  of  infinite  length  in  a  finite  time.  The  term 
fractal  was  invented  by  B.  Mandelbrot  (1977)  to  describe  a  class  of 
continuous  plane  contours  which  have  Infinite  length  and  no  derivatives; 
examples  of  (almost)  fractals  are  shown  in  Figs.  15  and  16  (taken  from 
Mandelbrot,  1977).  Given  some  plausible  assumptions  listed  below,  it  turns 
out  that  the  power  spectrum  of  a  plasma  cloud  whose  outer  contour  is  a 
fractal  is  reasonably  well  described  by  the  Rufenach  power  spectrum  (see 
equation  (II-7) 


P(k) 


Const. 


(I  + 


(III-3) 


with  L  (the  outer  scale)  of  the  order  of  the  overall  cloud  size  and  with 
g  significantly  le;  ■  than  2.  The  consequences  of  a  <  2  have  already 
been  pointed  out  (see  Section  II.  C). 
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15  Evolution  of  a  fractal  (in  thi*  case,  the  Koch  triadic) 
Taken  from  Mandelbrot  (1977). 
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For  a  geometric  fractal  like  those  in  Figs.  15  and  16,  a  is  a 
unique  number,  but  for  actual  plasma  clouds  there  will  be  an  ensemble  of 
a's  clustered  around  a  mean  value.  It  is  an  important  conclusion  our 
investigation  that  these  o's  — which  refer  to  a  fully-f ractalized  cloud — 
can  be  determined  from  computer  simulations  that  have  a  dynamic  range  no 
greater  than  now  available  {“  30)  . 

No  physical  cloud  can  ever  fully  fractalize,  because  at  a 
sufficiently  small  spatial  scale  there  is  physics  not  included  in  the 
canonical  equations  (e.g.,  electron  diffusion),  which  destroys  the  self¬ 
similar  evolution  of  the  cloud,  as  we  have  discussed  in  Sec.  Ill  A.  But  it 
is  worth  knowing  that  an  idealized  scale-free  cloud  can  fractalize  in  a 
finite  time,  of  the  order  of  twice  the  tim;.  it  takes  for  the  first 
generation  of  striations  to  appear.  This  is  because  each  new  strlation 
generation  evolves  on  a  time  scale  faster  than  that  of  the  previous 
generation,  in  proportion  to  the  size  ratio  between  the  old  and  the  new 
generation. 


1978). 


1 .  The  Canonical  Equations  and  Fractallzatlon 

These  are  (e.g.,  Perkins  et  al.  1973'  McDonald  et  al. 


3n  .  C( 


3t  f  •  ^n  »  0 


(III-4) 


^  •  (nf^)  =  ^  ^ 


(III-5) 
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Here  the  magnetic  field  B  points  in  the  z  direction,  4)  is  an 
effective  electrostatic  potential  Induced  by  the  plasma  cloud  of  (ion  or 
electron)  density  n  ,  and  is  an  externally-imposed  electric  field  in 

the  xy  plane.  The  cloud  potential  4>  vanishes  at  spatial  Infinity; 
summarizes  the  effects  of  ion-neutral  drag  and  of  gravity.  Ihe  fields  n  , 
({)  are  Independent  of  z  . 

Long  ago  Dungey  (1958)  proved  that  these  equations  have  no 
equilibrium  solutions  (aside  from  some  special  cases  which  will  not  concern 
us  here).  Again,  aside  from  some  special  cases,  there  are  no  known 
analytic  time-dependent  solutions,  but  It  is  well-understood  that  dense 
parts  of  the  cloud  tend  to  short  out  the  background  electric  field  and 

thus  move  more  slowly  than  the  background  drift  velocity  Uq  *  ^EqB  ^  .  As 
a  result,  the  back  side  of  the  cloud  tends  to  steepen,  then  to  split  into 
fingers.  The  NRL  simulations  show  that  when  the  first  few  backside  fingers 
have  formed,  the  edges  of  the  fingers  have  density  gradients  perhaps  an 
order  of  magnitude  steeper  than  the  Interior  of  the  fingers  and  of  the  rest 
of  the  cloud;  these  gradients  become  larger  as  time  goes  on.  In  view  of 
this,  it  has  been  suggested  that  the  power  spectrum  is  dominated  by  the 
steep  edges  (Zabusky  and  Block,  1978,  and  Costa  and  Kelley,  1978). 

For  a  steep-edged  plasma  cloud  of  well-behaved  shape  in  the  xy 

2  —2 

planes  and  typical  dimension  L  ,  this  means  P(k  )  ~  k  when  kL  »  1 
(but  k”^  is  greater  than  the  edge  thickness).  But  there  is  no  reason  to 
suppose  that  the  cloud  does  have  a  well-behaved  figure,  simply  because 


there  is  no  Intrinsic  scale  length.  A  roughly  circular  cloud  of  dimension 
L  will  develop  fingers  of  dimension  (1/2-1/3)  L  (say)  in  a  time  scaled 
by  T  =  (however,  the  time  scale  may  be  a  fraction  of  the  actual 

time).  If  there  are  no  collective  mechanisms  by  which  one  finger  interacts 
with  remote  parts  of  the  cloud  (and  none  has  been  demonstrated),  each 
finger  will  develop  fingers  just  as  the  original  cloud  did;  but  now  with  a 
time  scaled  by  (1/2-1/3)  LUo~^  ♦  It  is  easy  to  see  that  the  cloud 
develops  a  spectrum  of  finger  sizes  reaching  all  the  way  down  to  zero 
length,  and  it  does  this  in  a  finite  time,  measured  on  a  scale 
(1.5  -  2)LUq~^  .  Furthermore,  it  is  likely  that  in  this  idealized  cloud 
the  cloud  length  (measured  along  a  specific  isodensity  contour  at  the  edge) 
is  infinite  (of  course,  its  area  is  unchanged  since  the  velocity  field  is 
incompressible) . 

Mandelbrot  has  termed  such  geometric  oddities  fractals.  They  do 

not  have  the  same  dimension  as  ordinary  curves.  Here  the  dimension  is 

defined,  not  in  the  usual  topological  sense,  but  as  the  Hausdorff- 

Eesicovitch  dimension  (which  is  1  for  ordinary  curves).  Nor  is  the  power 

_2 

spectrum  of  a  steep-edged  cloud  whose  cross-section  is  a  fractal  »  k 
for  kL  »  1  ;  instead,  it  is  roughly  of  the  form  (I1I-3)  with  a  <  2  . 

(a  is  not,  unfortunately,  the  same  as  the  Hausdorff-Besicovitch 
dimension. ) 

In  the  rest  of  this  Section  we  1)  discuss  a  special  geometry  for 
equations  (III-4)  and  (III-5)  which  hints  at  fractallzatlon;  2)  coTiq)ute  the 
power  spectrum  of  some  fractals,  and  3)  deduce  from  2)  the  nhysical 


Ingredients  necessary  to  calculate  a  in  the  general  case.  These 
Ingredients  can  be  read  off  a  computer  simulation  with  limited  dynamical 
range. 

The  physics  is  cartooned  in  Fig.  17,  showing  how  a  plasma  cloud 
which  has  already  developed  steep  (backside)  edges  and  fingers  begins  to 
evolve  smaller  fingers  in  a  self-similar  way.  We  have  idealized  the  cloud 
to  have  a  uniform  density  except  for  an  infinitely  steep  edge  where  the 
density  drops  to  ambient;  this  is  a  popular  and  well- justified 
idealization,  judging  from  the  simulation  results.  Equation  (III-5)  shows 
that  polarization  charges  develop  on  the  sides,  but  not  the  backs,  of  the 
striations,  in  the  approximation  t  .  These  polarization  changes 

distort  Eq  by  an  amount  6^  ,  producing  wiggly  field  lines  (two  of  which 
are  shown  as  dotted  lines  in  Fig.  17).  The  velocity  flow  field 
5  =  cB~^Z  X  carries  plasma  out  of  the  center  of  the  striation, 

and  toward  the  edges.  As  a  result,  the  density-depleted  ceiter  cannot 
short  out  the  electric  field  as  well  as  the  edges  can,  and  the  center  moves 
more  rapidly  (to  the  right,  in  Fig.  17)  than  the  edges  do.  So  each  finger 
develops  an  indentation,  ultimately  resulting  in  two  new  fingers  of  about 
1/3  the  size  of  the  original  ones. 

The  idealized  outcome  of  several  such  generations  of  finger 
creation  is  shown  in  Fig.  18,  which  Illustrates  an  evolving  fractal.  This 
many  generations  would  be  calculable  in  standard  simulations  as  far  as 
spatial  dynamic  range  goes,  but  ordinarily  time  constraints  allow  at  most 
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Figure  17  Perturtwd  electric  field  ceueed  by,  erxf  cauting,  strietions. 

The  beckgrouFKl  pletme  it  to  die  left  and  doud  to  the  right. 
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two  generations  to  be  followed  [see,  for  example.  Figs.  2,  3,  4  of  Ossakow 
et  al.  (1977)]. 

By  now  it  is  evident  that  the  peculiar  structure  of  a  fractal 
looks  the  same  no  matter  what  power  magnifier  is  used  to  study  it. 
Everything  that  can  be  said  about  a  fractal  follows  from  knowledge  of  how 
one  generation  evolves  into  the  next,  which  requires  study  of  a  rather 
limited  dynamical  range.  We  show  how  this  works  in  determining  the  power 
spectrum  of  a  fractal. 

The  power  spectrum  of  equation  (111-3)  is  defined  as 

P(k)  -  ln(k)l2  (III-6) 


where 


n(k) 


f 

I  dx  dy  e 


Ik  ‘x 


(III-7) 


is  the  Fourier  transform  of  the  density  n  .  We  assume  n  is  constant 
inside  a  certain  isodenslty  contour,  dropping  to  (nearly)  zero  outside,  and 
take  ^  «  (k,0,,0)  (thus  Jc  •  5  «  0)  .  Consider  first  n(k)  for  a  non¬ 
fractal  shape,  such  as  the  topmost  triangle  in  Fig.  15.  Let  the  length  of 
a  side  be  L  .  For  kL  «  1,  n(k)  »  n(0)  *  A  ,  where  A  ~  L^  is  the  area 
of  the  figure,  and  for  kL  »  1  n(k)  ~  k~^  times  an  oscillating  function 
like  sin  kL  .  Then  P(k)  ~  k“^  ,  after  averaging  out  rapid  oscillations 
(sin^kL  +  1/2)  .  Then  for  a  non-fractal  shape  a  *  2  ,  which  we  call  the 


canonical  value  of  a  ;  it  is  often  invoked  as  a  rough  approximation  for 
actual  plasma  clouds  (Zabusky  and  Block,  1978,  and  Ck>sta  and  Kelley,  1978). 

Now  consider  the  power  spectrum  of  a  fractal,  for  a  given  kL  . 

In  Fig.  15  the  fractal  is  formed  by  adding  triangles  of  1/3  size  to  the 

existing  triangle  of  any  generation;  similarly,  n(k)  is  formed  additively 

•V 

from  triangles  of  size  L  x  3  ‘  ,  where  M  is  the  generation  number  (M  =  0 
for  the  original  triangle,  or  "unperturbed  cloud").  Let  us  choose  k  so 
that 

kL  =  3^  (III-8) 

for  some  Integer  M  rather  larger  than  1,  so  kL  »  1.  Then  in  computing 
the  contribution  to  n(k)  from  triangles  of  the  Mth  generation,  if 
M  »  M  the  exponent  ikx  in  (III-7)  hardly  varies  across  such  a 
triangle,  and  the  contribution  to  n(h)  is  just  Aj^  ,  the  area  of  the  Mth 
generation  triangles.  Similarly,  if  M  «  M  ,  ikx  varies  rapidly  and 
these  triangles  yield  a  term  in  n(k)  which  is  k~^  times  an  oscillating 
function.  We  drop  such  terms,  because  the  terms  we  save  do  not  decrease  as 
rapidly. 

It  is  a  reasonably  good  approximation  to  replace  the  inequality 
M  »  M  by  H  >  M  in  the  above  discussion,  which  leads  us  to  the 
approximation  for  kL  »  1  : 


MM  I 


OD 
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For  the  fractal  of  Fig.  15,  it  is  an  easy  calculation  to  find 


v-lv?) 


(M  >  1) 
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where  Aq  =  1/4  3L^  .  (Incidentally,  this  means  that  Fi®.  15  does  not 
quite  evolve  incoinpressibly ,  since  the  full  fractal  area  is  (8/5)  Aq»  but 
this  is  not  important;  it  is  quite  easy  to  change  Fig.  15  into  an 
incompressible  fractal  by  rescaling  the  whole  figure  a  slight  amount  at 
each  generation.)  Then  from  (III-9),  for  kL  »  1 


27  ^  -Fln(9/A) 
n(k)  .  Aj,e 
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and  using  (III-8), 


So  by  reference  to  the  Rufenach  spectrum  (711-3) 


(111-12) 


In  3 
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We  can  fit  the  power  spectrum  to  the  Rufenach  spectrum  at  another  point. 


say  k  =  0  ,  where  n(k)  «=  A  =  8/5  . 


This  shows  that  the  outer  scale  is 


(32/27)^^^^^  L  •>  h  ,  which  Is  intuitively  obvious:  the  sub-fractal 
generations  of  Fig.  15  do  not  change  their  overall  shape  very  much,  only 
their  edge  structure.  (This  last  point  is  related  to  the  phenomenon  of 
freezing;  it  shows  that  there  is  a  kind  of  shape  freezing  even  in  the 
presence  of  infinitely  detailed  bifurcation.) 

It  should  now  be  clear  how  to  calculate  a  for  any  fractal.  The 
formula  is 


a  -  2  ln|^^^j(lnX)~^  (III-14) 

where  is  the  area  of  the  Mth  generation  part  of  the  fractal,  and  X 

is  the  generation  ratio  (ratio  of  length  scales  of  the  Mth  to  the  (M  +  1) 
generation.  Applying  this  to  the  fractal  of  Fig.  18  (with  the  x  axis 
running  vertically)  gives  a  disturbingly  small  value  of  a  ,  about  0.75! 
This  is  because  the  generating  elements  in  Fig.  18  are  much  longer  than 
they  are  wide,  so  X  is  not  quite  well-defined.  Special  techniques 
involving  coupling  of  more  than  two  generations  are  needed,  but  we  will  not 
describe  them  here.  Nevertheless,  it  is  a  case  that  is  of  practical 
importance  judging  from  the  long  thin  fingers  of  the  con^uter  simulations. 

Since  no  physical  cloud  can  ever  fully  fractallze,  the  next 
question  is,  how  accurate  are  these  considerations?  One  can  redo  the 
problem,  replacing  the  upper  limit  M  in  (III-9)  by  a  finite  upper 
limit  My  ,  determined  by 
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L/£ 


(m-15) 


where  I  is  the  inner  scale,  or  scale  of  physical  effects  left  out  of  the 
canonical  equations  (e.'.,  £  «  Rl)  .  If  L/£  *  10^  ,  X  ■  3  ,  then 

•>  13  ,  and  the  error  committed  by  using  an  upper  limit  of  »  is  (for 

M 

the  Fig.  lb  fractal)  of  order  (^/9)  “  ~  3  x  10“^  as  long  as  kL  »  1  , 
k£  «  1  .  For  k£  >  1  ,  there  is  a  peak  in  the  power  spectrum  which  now 
begins  to  fall  at  least  as  k”^  and  possibly  faster,  depending  on  the 
details  of  how  steep  the  densit  gradient  is  on  the  scale  length  £  . 

These  questions  are  taken  up  in  Sec.  Ill  A. 


C.  Comparison  of  Instability  Evolution  in  Fluid  and  Particle 
Computer  Simulations  of  a  Plasma 

Until  rfew  }ntly  there  have  been,  in  effect,  two  different 
kinds  of  computer  simulations:  fluid  codes  (like  that  used  by  NRL)  which 
handle  complicated  spatial  structure,  but  not  particle-kinetic  effects;  and 
particle  codes,  used  primarily  for  infinite  homogeneous  plasmas.  Now, 
however,  it  is  possible  to  follow,  to  some  extent,  the  evolution  of 
gradients  in  particle-kinetic  codes,  as  John  Dawson  and  his  colleagues  at 
UCLA  have  shown. 

A  group  at  UCLA  under  the  direction  of  John  Dawson  has  done 
simulations  in  the  following  geometry  (see  Fig.  19,  taken  from  Goede's 
thesis,  Goede,  1980):  a  line  current  I  in  the  y  direction  creates  a 
magnetic  field  decreasing  as  x  increases.  A  plasma  cloud  (of 
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particle  "rods"  extending  along  fi  )  is  supplied,  also  with  a  gradient  in 

the  negative  x  direction;  of  course,  this  is  Rayleigh-Taylor  unstable. 

The  plasma  is  collisionless;  the  ion/electron  mass  ratio  M^/Mg  *  4  ,  and 

Tg/Tj^  =  2  .  The  evolution  of  the  particle  cloud  is  shown  at  four 

successive  times,  in  Fig.  20a,b,c,d.  The  cloud  moves  in  the  positive  x 

direction,  as  would  be  expected  from  the  linear  Rayleigh~Taylor  instability 

(with  "gravity"  g  given  by  -pBM  ^  ,  where  p  is  the  magnetic  moment), 

and  reaches  a  terminal  velocity  ~gv  ^  .  The  effective  collision 

eff 

frequency  is  not  yet  completely  understood,  but  it  may  arise  from  the 
generation  of  convective  cells  (Okuda  and  Dawson,  1973).  Computer 
printouts  of  the  electrostatic  potential  associated  with  Fig.  20  are  not 
yet  available,  but  printouts  for  similar  runs  in  Goede's  thesis  show  the 
formation  of  convective  cells  which  appear  to  be  responsible  for  the 
Induced  drag  on  the  moving  blob  of  plasma. 

The  s^Lraulations  show  unstable  electrostatic  modes,  but  the  growth 
rates  of  these  modes  is  not  at  all  well-fit  by  simple  collisionless 
theory.  We  have  pointed  out  to  Goede  and  Dawson  that  the  growth  rates  may 
possibly  be  fit  by  a  collisional  Rayleigh-Taylor  form,  using 
given  from  the  terminal  blob  velocity. 

For  those  interested  in  striations,  the  most  interesting  aspect  of 
the  simulations  is  the  formation  of  fingers  on  the  blob  (Fig.  2&d).  These 
fingers  are  precisely  of  the  type  expected  to  form  as  a  result  of  the 
2  X  S-gradient  drift  instability,  which  is  mathematically  the  same  as  the 
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One  of  four  snapshots  of  a  plasma  cloud 
of  Fig.  19.  Time  Is  in  units  of  the  el 
and  M./M  “A,  (Furnished  by  .lohn  Daws 


collislonal  Rayleigh-Taylor  Instability  (with  effective  electric  field 
E  •=  (M/e)g(n/Vg^p)  . 

The  way  now  lies  open  to  combine  large-scale  fluid  codes  and 
smaller-scale  particle  codes,  since  there  is  a  real  overlap  in  their 
domains  of  applicability.  Dawson  and  colleagues'  results  Indicate  that 
convective  cells  are  formed;  in  the  presently  available  simulations,  cell 
formation  is  enhanced  by  having  a  strong  magnetic  field 
will  be  most  interesting  to  stud”  weak-field  simulations,  looking  for 
convective  cell  generation  caused  by  the  spontaneous  formation  of  fingers 
(i.e.,  relatively  sharp  gradients).  We  have  already  indicated  in 
Sec.  Ill  A  how  this  might  work,  through  the  mechanism  of  drift-wave 
generation  and  decay  into  convective  cells.  The  proposed  analytical 
research  of  such  a  mechanism  must  be  accompanied  by  simulations  capable 
simultaneously  of  dealing  with  particle-kinetic  effects  and  with 
spontaneous  striatlon  formation  by  the  standard  electrostatic-fluid 
Instabilities;  these  simulations  now  appear  to  be  possible. 

D.  Possible  Hagnetospheric  Fffect*^ 

A  high-altitude  nuclear  explosion  generates  an  enormous 
number  of  energetic  ions.  It  is  conventional  wisdom  that  these  are 
eventually  cooled  by  collisions,  and  fall  back  into  the  ionosphere. 
However,  the  possibility  exists  that  substantial  numbers  of  energetic  (>  1 
kev)  ions  undergo  pitch-angle  scattering  at  very  high  altitudes  which  puts 
them  into  orbits  whose  mirror  points  are  so  high  that  they  are  virtually 
collisionless.  (The  pitch-angle  scattering  may  be  caused  by,  e.g.. 
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electromagnetic  ion-cyclotron  waves  generated  by  the  energetic  ions 
themselves.)  Ihese  ions  may  live  for  tens  of  hours  or  days,  their 
lifetimes  being  governed  by  charge  exchange,  pitch-angle  scattering  into 
the  loss  cone,  and  very  slow  energy  degradation  by  Coulomb  collisions. 


Such  energetic  ions  are  impossible  to  detect  directly  from  the 
ground  (by,  e.g.,  optical  observation).  As  far  as  we  know,  there  are  no 
satellite  data  which  bear  one  way  or  the  other  on  trapped  Ions  from  a  high- 
altitude  nuclear  explosion,  because  of  lack  of  adequate  Instrumentation  in 
the  years  1963  and  before;  nor  is  there  any  theoretical  work  (including 
simulations)  on  the  problem,  the  direct  scintillation  effects  of  the  ion 
density  are  clearly  negligible  because  the  densities  are  so  small.  Here  we 
are  looking  for  effects  in  the  more  dense,  lower  ionosphere  which  are 
caused  by  the  interaction  of  these  ion  currents  with  the  ionospheric 
plasma. 
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If  such  a  partial  ring  current  of  energetic  ions  is  found,  it  is 
worth  looking  at  its  effect  on  the  striation  scene.  Ihis  ring  current  will 
Itself  have  sharp  gradients  which  could  influence  propagation  from 
satellites.  Moreover,  these  sharp  gradients  will  drive  field-aligned 
currents  which  can  Interact  with  conventional  ionospheric  striation 
mechanisms  (e.g.,  by  driving  a  current-convective  instability,  or  by 
generating  high-altitude  electrostatic  turbulence  which  is  mapped  down 
field  lines  into  the  ionosphere).  Ihese  currents  come  from  gradient  and 
curvature  drift  velocities  of  some  hundreds  of  meters  to  a  few  km/second. 
One  might  speculate  that  the  artificial  ring  current  will  generate  a  kind 
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of  auroral  environment,  associated  with  (small)  electric  fields  parallel  to 
the  magnetic  field.  (This  speculation  rests  on  recent  auroral  models  which 
drive  parallel  electric  fields  by  energetic  ion  Injection.)  The  resulting 
currents  and  energetic  electrons  may  influence  striation  phenomena  by 
providing  rapid  cross-field  electron  transport,  thus  allowing  Ions  to 
diffuse  (see  Sec.  Ill  A).  It  is  worth  a  modest  theoretical  effort  to  see 
if  any  such  ion  ring  current  is  actually  formed,  and  if  so,  what  its 
effects  on  strlatlons  are. 

It  has  recently  been  established  in  an  experiment  called  Trigger 
that  auroral  phenomena,  including  rapid  acceleration  of  particles  to 
energies  of  as  much  as  300  keV,  can  be  stimulated  by  the  s»>'’den  release  of 
large  amounts  of  plasma  in  the  F  region  ionosphere  (=>  160  km).  The 
Trigger  experiment  (Holmgren  et  al.,  1980)  consisted  of  an  explosive 
release  of  Cs"*"  io.is  and  electrons  (via  an  explosion  of  12  kg  of  a  TNT-AIO- 
CSNO3  mixture)  which  increased  the  conductivity  by  more  than  2  orders  of 
magnitude  over  ambient.  During  the  few  seconds  to  few  tens  of  seconds  of 
observation,  the  cloud  expanded  to  a  scale  size  of  a  few  km.  There  were  no 
energetic  particles  to  begin  with  in  the  plasma;  apparently  they  were 
created  by  a  large  transient  electric  pulse  with  a  substantial  component 
along  the  magnetic  field  line.  Obviously  such  an  electric  field  will 
accelerate  ambient  ions  (a  process  which  is  not  necessarily  modeled 
correctly  in  MHD  codes,  since  kinetic  effects  may  dominate  the  parallel 
electric  fields),  so  that  it  may  well  be  that  energetic  ions  in  the 
magnetosphere  are  not  simply  bomb  debris,  but  strongly-energized 
ionospheric  ions. 


One  may  well  ask:  Should  not  such  phenomena  have  been  seen  In  Ba*^ 
releases?  In  fact,  it  would  not  be  surprising  to  see  them,  and  there  is 
some  evidence  that  they  were  (see  Baker  et  al.,  1978,  for  an  overview),  but 
until  Trigger  little  attention  was  paid  to  proper  instrumentation  and 
especially  to  Insuring  that  instruments  were  on  the  same  field  line  as  the 
released  plasma. 

It  would  be  very  much  worthwhile  to  conduct  Trigger-like 
experiments  in  future  Ba"*"  releases,  not  only  for  what  is  learned  about 
early-time  processes,  but  for  what  might  be  learned  about  the  storage  and 
subsequent  release  at  late  times  of  energetic  magnetospherically-trapped 
particles. 
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IV.  LABORATORY  EXPERIMENTS  AS  MODELS  OF 

HIGH  ALTITUDE  NUCLEAR  EXPLOSION  PHYSICS 

In  this  section,  the  value  of  laboratory  experiments,  specifically 
laser  exploded  pellet  experiments,  to  obtaining  a  more  complete 
understanding  of  the  phenomenology  of  a  high  altitude  nuclear  explosion 
(HANE)  is  discussed.  For  the  reasons  discussed  below,  it  appears  that  such 
experiments  can  provide  Information  regarding  (1)  the  "early  time"  coupling 
of  the  expanding  debris  cloud  to  the  ambient  ionosphere;  (2)  the  importance 
of  "early  time"  effects  in  determining  the  later  time  phenomena;  and  (3) 
the  late  time  inner  and  outer  scale  sizes.  This  information,  and  the 
u''de-  standing  of  HANE  phenomenology  to  which  it  would  contribute,  would 
come  both  from  the  results  of  the  experiments  themselves  and  from  the 
improvement  that  would  undoubtedly  follow  in  the  computer  codes  used  to 
simulate  HANE  physics.  The  benefits  of  improved  ability  to  predict  the 
long  term  behavior  of  the  ionosphere  after  a  HANE  could  have  a  substantial 
impact  on  the  effectiveness  of  our  satellite  communication  system,  as 
discussed  in  detail  in  Section  II. 

In  the  subsections  which  follow,  we  further  motivate  HANE 
simulation  with  laser  pellet  explosions,  briefly  describe  experiments  which 
have  been  performed  or  are  in  progress,  and  then  discuss  in  some  detail  the 
possibility  that  relevant  experiments  can  be  performed  using  existing 
facilities  with  relatively  modest  additional  capital  cost.  We  shall  see 
that  two  independent  experiments  are  probably  necessary  to  simulate  early 
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and  late  cine  phenomenology.  It  is  important  to  note  chat  we  are  not 
providing  an  experimental  design,  but  merely  a  strong  Indication  that 
relevant  experiments  can  be  designed.  In  order  to  be  sure  that  proper 
consideration  of  all  necessary  scaling  laws  is  given,  it  is  recommended 
that  DNA  support  two  to  three  month  design  studies  by  interested 
organizations  who  have  the  necessary  major  pieces  of  equipment  (powerful 
lasers,  sophisticated  diagnostic  capability,  and/or  large  vacuum 
chambers).  This  would  enable  DNA  to  evaluate  whether  truly  relevant 
exTieriments  can  ba  performed,  after  which  an  RFP  for  such  experiments  car. 
be  Issued. 

A.  Objectives  of  the  Suggested  Experiments 

The  modeling  of  HANE  phenomena  with  barium  releases  in  the 
ionosphere,  together  with  analytic  theory  and  computer  simulation,  has  led 
to  considerable  understanding  of  the  fundamental  physics  underlying  the 
behavior  of  artificial  plasma  clouds  in  the  ionosphere  during  the 
electrostatic  regime  (see  Section  I).  (Evolution  at  very  late  times  is 
still  not  understood  either  empirically  or  theoretically;  see  Sections  I 
and  Til.)  Nevertheless,  our  inability  to  make  continuous  detailed 
measurements  on  barium  clouds,  and  the  fact  that  barium  clouds  can  simulate 
only  electrostatic  regime  effects,  leave  us  uncertain  as  to  (1)  the 
evolution  of  a  plasma  cloud  in  the  ionosphere,  specifically  its  inner  and 
outer  scale  lengths,  and  the  spatial  power  spectrum  of  irregularities,  and 
(2)  the  degree  to  which  the  details  of  early  time  phenomenology  affect  the 
lace  time  behavior  of  HANE  plasma  clouds.  Furthermore,  we  are  unable  to 
confidently  infer  all  of  the  necessary  details  from  computer  simulation 
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coupled  with  theory.  This  is  true  even  in  time  regimes  where  we  have  codes 
which  are  "validated"  against  barium  cloud  experiments.  This  is  because 
limitations  on  the  number  of  spatial  grid  points  that  can  be  used  in  a 
simulation  of  adequate  time  duration  limit  the  range  of  spatial  scales  that 
can  be  covered  to  a  factor  of  about  100,  e.g.,  0.8-80  km  (see,  for  example, 
Kesklnen,  1980).  By  contrast,  in  an  actual  barium  release,  wavelengths 
ranging  from  <  10  m  to  about  50  km  appear  to  be  Important  (Baker  et  al. 
1978).  As  for  early  times,  there  is  still  controversy  over  the  actual 
debris-air  coupling  mechanlsm(s)  (Clark,  1980),  and  these  may  provide 
initial  conditions  for  subsequent  late  time  plasma  cloud  expansion  which 
substantially  Influence  the  Important  scale  sizes  and  the  power  spectrum  of 
the  late  time  cloud  structure. 

Given  the  above  background  it  is  clear  that  appropriate 
experiments  would  be  quite  valuable.  The  goals  would  be  the  following: 

(1)  an  experiment  which  investigates  the  early  time  debris-air  coupling  and 
would  thus  provide  information  that  can  be  otherwise  provided  only  by  a 
HANE;  this  experiment  would  also  serve  to  provide  information  on  the 
initial  conditions  which  should  be  used  in  late  time  computer  simulations; 

(2)  an  experiment  which  investigates  the  late  time  plasma  cloud  behavior, 
and  which  can  be  used  to  determine  relevant  length  scales  and  to  "validate" 
late  time  codes  if  it  is  well  diagnosed*  As  we  will  see,  it  would  be 
wasteful  to  tie  up  a  laser  capable  of  doing  the  first  experiment  to  also  do 
the  second  one.  Therefore,  two  independent  experiments  are  suggested 
below,  one  using  a  very  high  power  laser  and  the  other  using  a  lower  power 
laser,  but  a  very  large  vacuum  chamber. 


B.  Closely  Related  Laboratory  Experiments 

About  ten  years  ago  a  laser-pellet  experiment  for  HANE 
simulation  was  initiated  at  the  Naval  Research  Laboratory  (NRL)  by  John 
Stamper  and  Steve  Dean.  A  10-60J,  30  nsec  pulse  duration  Nd:glass  laser 
was  focussed  on  the  tip  of  a  250  urn  diameter  Incite  fibre.  A  1  kG  magnetic 
field  was  to  be  used,  as  was  a  10^^  cm“^  "ambient"  plasma  (Stamper, 

1980).  However,  the  background  plasma  was  found  to  be  unnecessary  because 
the  laser-produced  plasma  at  the  target  was  able  to  substantially  photo- 
lonlze  10-200  mtorr  of  nitrogen  or  helium  in  the  vicinity  of  the  target. 

The  intent  was  to  look  at  the  early  time  "collisionless"  portion  of  the 
interaction.  The  experiment  was  sidetracked  by  the  discovery  and  study  of 
large  spontaneous  magnetic  fields  near  the  pellet  (Stamper  et  al»  1971). 
These  fields,  which  were  due  to  the  single-sided  pellet  Irradiation  at  high 
power  density  (-  10^^  W/cm^),  dominated  the  early  time  plasma  expansion. 

At  the  radius  at  which  the  spontaneous  fields  were  no  longer  dominant,  the 
electrons  were  no  longer  collisionless,  although  the  ions  were  (Dean  et  al. 
1971).  No  follow-on  experiment  with  symmetric  irradiation  of  the  pellet 
was  performed  at  NRL. 

At  the  University  of  Maryland,  an  experiment  to  follow  up  the  work 
of  Stamper  and  Dean  on  HANE  simulation  was  started  in  1971.  The  problem  of 
the  megagauss  fields  was  avoided  by  lowering  the  laser  power  density, 
looking  at  least  1  cm  away  from  the  target,  and  using  a  lower  ambient 
plasma  density.  Early  experiments  found  a  collisional  interaction  (Koopman 
and  Goforth,  1974),  but  in  more  recent  experiments,  the  interaction  Is 
believed  to  be  collisionless  (Parsons,  1980).  In  these  recent  experiments 
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a  slab  target  of  Ba  or  Iron  is  used,  and  it  is  being  irradiated  by  a  6J  CO2 
laser  in  ~  100  ns.  This  means  that  the  Initial  expansion  velocity  is  much 
lower  than  in  a  HANE,  namely  an  observed  value  of  (1-3)  x  10^  cm/sec. 
Observations  are  made  where  the  expanding  plasma  density  is  from  10^^  to 
10^3 

cm“^.  The  ambient  plasma  of  10^ ^  cm"^  is  produced  from  the  10"^ 

Torr  (6  x  10^^  cm~^)  background  gas  by  the  U.V.  flash  of  the  target.  The 
magnetic  field  is  variable  from  0-2  kG.  This  experiment  has  observed 
collisionless  coupling  with  1  cm  shock  thickness  using  Ba,  with  a  1  kG 
magnetic  field.  The  specific  instabilities  involved  in  the  collisionless 
interaction  have  not  yet  been  identified,  but  this  is  precisely  what  this 
experiment  is  now  capable  of  doing.  Apparently  the  experiment  has  been 
funded  at  a  level  of  ~  $80k/year  by  the  Air  Force,  during  which  its 
diagnostic  capability  was  built  up  to  a  very  formidable  level.  It  is  now  a 
beautiful  physics  experiment  capable  of  determining  which  of  the  standard 
collection  of  instabilities  are  present  when  one  plasma  expands  into 
another.  However,  the  experiment  has  been  terminated  due  to  the  recent 
death  of  the  principal  Investigator,  David  Koopman.  An  SAI  group  has 
proposed  to  AFOSR  that  the  experiment  be  continued  under  their  auspices. 

The  results  of  this  experiment  are  significant  to  the  physics  of 
the  HANE  problem.  The  findings  regarding  collisionless  coupling  in  the 
most  recent  experiments  should  be  examined  closely  for  potential 
applicability  and  Importance  to  a  HANE.  However,  the  specific  energy 
deposition  (J/gm)  is  much  lower  than  in  a  HANE  [Zinn  et  al.  (1965)  give  the 
kinetic  energy  yield  for  Starfish  of  about  lO^J/gm) ,  and  a  slab  target  was 
used  rather  than  a  pellet.  Since  specific  energy  deposition  is  said  to  be 
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one  of  the  most  important  scaling  parameters  (Clark ,  1980) ,  this  may  make 
any  conclusions  from  the  University  of  Maryland  experiment  difficult  to 
extrapolate  with  confidence  to  the  early  and  intermediate  time  regimes  of 
the  high  altitude  case.  The  recent  work  on  this  experiment  was  apparently 
under  the  direction  of  Dr.  Charles  Parsons,  and  the  experiment  was  being 
carried  out  by  G.  P.  Jelllson,  a  graduate  student. 

A  variety  of  other  slightly  related  experiments  have  been 
performed  which  we  mention  only  briefly:  two  experiments  to  fill  mirror 
machines  with  laser-produced  plasma  have  been  performed,  one  by  Alan  Haught 
at  United  Aircraft  and  one  by  John  Osher  at  LLL  ("Baseball"  experiment). 

An  experiment  by  Pechacek  and  Grelg  at  MRL  produces  a  plasma  in  a  cusp 
magnetic  field  using  a  laser-exploded  pellet.  None  of  these  experiments 
involved  an  ambient  plasma,  so  there  was  no  coupling  to  be  observed. 
However,  data  on  plasma  production  and  expansion  is  presumably  available. 

An  experiment  at  AVCO  to  Investigate  the  coupling  of  a  gun  plasma  to  an 
ambient  plasma  was  performed  about  ten  years  ago  by  Friedman  and  Patrick 
(1971).  This  experiment  demonstrated  collisionless  coupling  of  a  streaming 
plasma  to  a  magnetized  background  plasma.  The  streaming  direction, 
however,  was  strictly  parallel  to  an  applied  aagner.lc  field,  so  the 
relevance  to  a  HANE  is  marginal.  Nevertheless,  it  Is  worth  noting  that 
momentum  transfer  lengths  from  streaming  to  background  plasma  a  factor  of 
100  less  than  the  ion  Larmor  radius  were  observed. 


C.  Present  Status  of  Pulsed  Laser  Systems 


Multikilojoule  laser  sytems  exist  in  several  laboratories  in 


the  U.S.  now,  with  the  largest  development  programs  being  those  at  LASL 
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(CO2  laser  system)  and  LLL  (Nd:Glass  laser  system).  We  shall  use  the  SHIVA 
laser  system  at  LLL  as  the  basis  for  discussion  because  it  delivers  the 
most  energy  to  a  target,  namely  10  kJ  to  a  <  200  un  diameter  pellet.  In  a 
short  pulse  (<  100  ps),  about  25X  of  the  energy  Is  absorbed,  and  virtually 
100%  can  be  absorbed  in  longer  (>  200  ps)  pulses  (Manes  et  al.  1980; 
Lawrence  Livermore  Laboratories,  1979).  Furthermore,  a  pellet  can  be 
Irradiated  with  nearly  arbitrary  symmetry.  A  typical  exploding  pusher 
target  is  a  glass  shell  150  pm  in  diameter,  1  im  wall  thickness,  filled 
with  10  atmospheres  of  D-T  gas.  It  weighs  100  ngm,  and  it  is  irradiated  in 
the  short  pulse  mode  with  close  to  spherical  symmetry.  A  specific  energy 
deposition  of  over  10^  J/gm  is  implied,  and  6  x  10^  J/gm  is  documented.  We 
specifically  mention  exploding  pusher  targets  because  they  are  rather  well 
characterized  by  now.  For  example  (Manes  et  al.  1980),  the  expansion 
starts  at  (1-10)  x  10^  cm/sec,  depending  upon  the  power  on  the  silicon 
microballoon  target.  Furthermore,  “other  observables",  e.g.,  shell 
temperature  and  core  information,  are  well  predicted  by  code  calculations 
at  the  early  times  which  Interest  the  LLL  group. 

In  addition  to  the  laser  Itself  the  SHIVA  system  includes  a 
multiport,  80  cm  radius  experiment  chamber  and  a  very  sophisticated 
collection  of  diagnostics.  The  latter  is  particularly  suited  for  very  high 
density  and  temperature  observations,  but  some  of  it  could  be  useful  to 
investigate  coupling  phenomena  at  lower  density  if  an  ambient  plasma  weze 
present  in  the  chamber.  Target  fabrication  has  also  become  rather 
sophisticated  in  recent  years:  thin  shell,  raultlshell  and  multilayered 
targets  have  all  been  made  in  many  different  materials.  Therefore,  it  is 
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likely  that  appropriate  adjustment  of  laser  power  and  pellet  surface 
material  would  allow  the  pellet  to  produce  the  ambient  plasma  around  It  by 
photolonlzatlon  as  a  HANE  does,  while  still  varying  other  observables,  such 
as  the  Initial  expansion  velocity. 

The  great  advantage  of  a  laboratory  experiment  Is  evident  when  one 
realizes  that  the  SHIVA  system,  for  all  Its  complication,  fs  capable  of  two 
laser-pellet  experiments  at  full  power  and  with  full  diagnostic  complement 
per  day  (Manes  et  al.  1980).  Thus,  once  appropriate  parameters  have  been 
determined  and  set  up  experimentally,  and  the  diagnostics  are  set  up. 
Information  can  be  collected  (and  analyzed)  very  rapidly.  This  Information 
can  be  obtained  using  a  variety  of  different  pellets  and  laser  powers, 
depending  upon  the  effect  under  Investigation.  A  variety  of  different 
symmetries  can  also  be  Investigated. 

D.  Scaling  from  HANEs  to  Laser  Pellets 

The  scaling  laws  to  be  used  to  guide  the  design  of 
experiments  to  simulate  a  HANE  are  relatively  straight-forward  for  a  late 
time  experiment,  but  not  nearly  so  obvious  for  the  early  time.  There  are 
two  dominant  problems  In  the  early  time  case:  (1)  It  Is  not  possible  to 
obtain  consistent  scaling  of  all  parameters,  and  so  It  Is  necessary  to 
choose  a  set  which  makes  the  experiment  "as  relevant  as  possible";  and  (2) 
there  is  no  agreement  as  to  what  coupling  mechanisms  are  the  most  Important 
under  given  conditions  (e.g.,  altitude  and  yield),  and  different  scaling 
may  be  required  depending  upon  the  mechanlsias  expected.  (Here  Is  where  the 
results  of  a  physics  experiment,  such  as  the  University  of  Maryland  one 
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described  in  the  previous  subsection,  would  be  very  helpful.)  Since  the 
least  information  is  available  for  the  early  time,  it  Is  probably  the  most 
important  case  to  consider  here,  and  so  we  shall  look  at  It  first.  This 
case  will  also  be  able  to  provide  information  on  the  transition  to  late 
time,  the  Importance  of  which  has  already  been  discussed. 

1.  Early  Time  Scaling 

Table  II  presents  the  parameters  of  the  Starfish  Device 
exploded  at  400  km  altitude,  as  presented  by  Zinn  et  al.  (1965),  in  the 
column  marked  HANE.  Although  we  consider  Starfish  here,  we  recognize  that 
other  event  paramete.'s  are  probably  of  as  much  or  more  Interest  and  that 
these  would  be  considered  in  further  work.  Note  that  near  the  explosion 
the  ambient  air  Is  likely  to  be  fully  ionized  by  the  flash,  but  further 
away  the  plasma  density  will  be  close  to  ambient.  Thus,  two  different  sets 
of  plasma  frequencies  are  given.  Because  the  atmospheric  density  scale 
length  at  400  km  is  ~  60  km,  the  debris  expansion  and  shock  wave  are  not 
symmetric  about  the  burst  point  in  the  vertical  direction  even  before  the 
magnetic  field  has  much  of  an  effect.  Furthermore,  the  symmetry  is  even 
worse  when  the  debris  energy  has  been  degraded  sufficiently  that  the 
magnetic  field  begins  to  become  relevant. 

With  these  considerations  and  others  implied  by  the  parameters  in 
Table  II  in  mind,  we  specify  that  the  experiment  should  have  the  capability 
of  depositing  at  least  lOOJ  in  a  100  ngm  target  with  some  appropriate  level 
of  symmetry.  Provisions  to  puff  in  gas  to  produce  a  neutral  particle 
density  gradient,  and  to  produce  a  magnetic  field  oriented  appropriately  to 
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TABLE  II 


EARLY  TIME  PARAMETERS 


HAME 

Experiment 

Kinetic  Energy  Yield 

10l5  j 

100  J 

Mass 

10®  gn 

100  ngm 

Bubble  Radius 

~  200  km 

5  cm 

"Mass  Radius"* 

100  km  (~  10®  cm”3) 

2  cm  (~  10 cm-3 

Magnetic  Field 

300  G 

Initial  ion  Larmor  radius 

(lOkeV/Nucleon) 

29/Z  Km(Fe) 

30  A  cm 

Background 

Air  (10®  cm“^) 

Air  (>  10^^  cm”^) 

(1) 

ce 

10^  sec“^ 

2  X  10^0  sec"^ 

<*>  j 

cl 

lO(Fe)  sec“^ 

lO^f 

A 

V  (leV) 
e 

3  X  10^  sec”^ 

3  X  10^  sec"^ 

''in 

~  1  sec"^ 

10^  sec“^ 

0)  (Fully  Ionized 

pe 

10®  cm®  Ambient**) 

5  X  10®  sec“^ 

5  X  10^^  sec"^ 

(Fully  Ionized  Ambient) 

6  X  10^  sec”^ 

6  X  10®  sec"^ 

to  (Undisturbed 
pe 

Ambient  2  x  10^  cm~^) 

2  X  10^  sec“^ 

2  X  lOl®  sec"l 

lOpi  (Undisturbed  Ambient) 

2  X  103  sec"^ 

2  X  10^  sec'^ 

(Undisturbed  Ambient) 

~  400  km/sec 

~  400  km/sec 

(Fully  Ionized  Ambient) 

~  20  km/sec 

~  20  km/sec 

*  The  radius  within  which  there  Is  an  ambient  mass  equal  to  the  weapon 
mass. 

**  Assuming  the  background  gas  is  fully  Ionized  by  the  radiation  in  the 
region  of  Interest* 
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that  density  gradient,  should  be  made,  although  the  first  experiments  could 

be  performed  without  these.  A  magnetic  field  of  up  to  ~  1  kG  is  necessary 

to  reduce  length  scales  to  laboratory  size.  This  is  seen  by  assuming  5^  x 

Volume  is  proportional  to  the  kinetic  energy  yield  (i.e.,  the  bubble  size 

is  determined  oy  field  energy  displacement  by  the  expanding  debris).  Thus, 

if  we  take  the  desired  laboratory  length  to  be  5  cm  and  the  remaining 

quantities  from  Table  II,  we  obtain  300  G  as  our  field.  In  order  to  set 

the  density,  we  have  two  possible  scaling  choices:  B^/n  or  B/n  the  same. 

In  addition  we  wish  to  have  the  "mass  radius"  approximately  equal  to  }/2 

the  bubble  radius.  Ihe  first  density  scaling  retains  the  value  of 

u^/Wpe  »  the  Alftfin  speed  Vy^  (assuming  the  same  ion  species)  and 

transverse  pressure  balance  (at  the  rame  temperature);  it  requires  a  plasma 

density  of  10^^  cm~^.  The  B/n  scaling  maintains  u  /v  (at  the  same 

ce  e 

temperature)  as  well  as  the  ratio  of  the  ion  Larmor  radius  to  the  ion 
range;  the  density  would  be  only  3  x  10^^  cm~^  here.  The  mass  radius  would 
be  2  cm  for  10^^  cm“^  air,  requiring  a  fill  pressure  of  1.5  mtorr  of  air  if 
the  U.V.  flash  fully  (singly)  ionizes  it.  We  choose  10^^  cm“^,  although  it 
may  be  necessary  to  drop  it  by  an  order  of  magnitude  to  avoid 
collisionality  problems —  table  values  may  be  too  large.  The 

other  parameter  poorly  scaled  in  the  experiment  is  the  ion  Larmor  radius, 
which  may  cause  a  problem  in  trying  to  use  the  results  of  this  experiment 
to  provide  initial  conditions  for  a  late  time  code  calculation.  This 
presumably  depends  upon  the  coupling  mechanism  and  requires  further 
study.  A  small  amount  of  prelonlzation  can  be  used  to  make  an  undisturbed 
ambient  density  of  2  x  cm~^  in  order  to  have  undisturbed  ambient  Vy^, 

and  u  scale  correctly.  The  interesting  duration  of  this  experiment 


IV-11 


v7ould  be  ~  10  (jsec  at  most*  Extensive  optical  diagnostics  will  probably  be 
needed,  such  as  holographic  Interferometry  for  density  and  s,>ectroscopy  for 
temperature*  Submllllneter  scale  measurements  can  be  made* 


(a)  T  *  experiment  duration  >  2  x  10^  “cl~^  order  to 
simulate  ~  1000  sec  Is  the  Ion  cyclotron  frequency) ; 

(b)  T  >  100  where  y  *  growth  rate  of  the  gradient 
drift  Instability  (c£/BL,  where  L  »  density  scale  length,  E  and  B  are  the 
electric  and  magnetic  fields,  respectively,  and  c  Is  the  speed  of  light); 


We  assume  we  have  a  >  S  m  diameter  by  many  meter  long  tank 
available  In  which  to  do  the  experiment.  Then  the  following  set  of 
experimental  conditions  will  suffice: 


I 

p 


t.:* 

it 


‘■j 


<  jS 


“2?  * 


i  ^ 


B  -  1  kG 
Ti  -  1000°  K 

E  ■=  50  V/n  across  the  5  lu  tank 


These  give  =  1.6  x  10^  sec~^,  ion  Larmor  radius  ~  1  mm,  and  a  time  of 
drift  across  the  5m  tank  of  10  msec.  In  order  to  have 
y  *  cE/BL  >  10^  sec“^  requires  L  <  5  cm.  The  densities  must  be  determined 
by  the  requirements  on  the  collision  and  plasma  frequencies.  In  addition, 
we  must  be  sure  that  the  Li"*"^  plasma  is  cooled  by  collisions  with  a  cold 
neutral  background  in  ^  100  jisec  (<  .  The  following  set  of 

parameters  seems  reasonable  based  upon  what  the  University  of  Maryland 
experiment  is  doing  with  laser  produced  barium  plasmas  from  a  slab  target; 
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n  >  6  X  10  cm  (>  1  mTorr) 
n  ~  " 


12  -3 

n„  «  10  cm 


10  -3  11  -3 

n  ,  «  10  cm  or  10  cm 
emb 


The  two  choices  in  the  ambient  density  are  to  simulate  BANE  and  barium 
cloud  situations,  respectively.  Assuming  the  ambient  plasma  has  cold  ions 
(<  0.1  eV)  and  electrons  (<  i  eV)»  all  of  the  necessary  frequency  criteria 
are  satisfied. 

Such  a  situation  may  be  produced  by  a  modest  power  laser.  The 
necessary  lO^cm^  of  10^2  li+I  is  lo^^e-i  pairs.  At  10  eV  to  produce 
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each  one,  this  requires  160  mJ.  A  slab  target  ~  1  cn  in  radius  vrould  be 


used.  The  plasma  would  expand  into  the  weakly  ionized  ambient  plasma  with 


an  initial  ion  Larraor  radius  of  ~  1  cm  (10  eV)  and  a  collision  frequency  of 


>  10^(10  eV).  Therefore,  in  100  psec  it  will  have  collided  many  times  and 


diffused  across  B  u  few  ion  Larmor  radii,  it  is  now  low  B  and  "ready* 


to  respond  to  a  250  V  electric  field  across  the  tank  perpendicular  to  3. 


The  tank  must  be  much  more  than  5  m  long  so  that  plasma  does  not  reach  the 


ends  during  the  10  msec  experiment.  Such  a  tank  exists  at  the  Manned 


Spacecraft  Center  near  Houston  (it  was  used  for  Space  Shuttle  system  tests) 


and  others  may  exist  elsewhere  (e.g.,  one  of  the  "space  chambers"  that 


exist  around  the  country).  (A  great  deal  of  capital  expense  can  be  saved 


if  a  tank  and  a  field  coil  power  supply  are  available  at  the  same  site. 


The  coils  themselves  and  the  laser  would  be  relatively  inexpensive.)  Note 


that  it  may  be  possible  to  introduce  perturbations  into  the  Initial  cloud 


by  the  use  of  periodic  barriers  in  the  way  of  the  plasma  as  it  moves  away 


from  the  slab,  or  perhaps  by  using  several  separated  focal  spots  for  the 


laser.  Since  assembling  the  necessary  diagnostic  capability  is  probably 


the  most  time-consuming  aspect  of  building  up  this  experiment,  a  reasonable 


approach  is  to  transport  the  University  of  Maryland  experiment  to  wherever 


an  adequately  large  tank  is  available.  (Even  if  this  is  not  done,  this 


experiment  should  be  kept  '^oing  in  order  to  complete  its  collisionless 


coupling  physics  experiments.) 


E.  Conclusion  Regarding  Laboratory  Experiments 


The  discussion  of  the  previous  section  by  no  means 


constitutes  an  experimental  design.  However,  it  does  indicate  that  a 
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complete  design  study  Is  warranted  if  DM  is  Interested  in  using  controlled 
and  flexible  laboratory  experiments  to  simulate  high  altitude  nuclear 
explosions.  Recommendations  are  made  In  Section  I  above. 


V.  RADIATIVE  TRANSFER  EFFECTS  IN  Ba+  LINES  IN 
STRIATED  BARIUM  CLOUDS 

2  2 

Ba  II  has  two  resonance  transitions,  6  ^1/2  and 

2  2 

6  Sj^,-6  ^212  slightly  longer  wavelengths  lie  three 

fluorescent  transitions,  5  ^^1/2  (^^97A),  5  ^^3/2 

2  2 

(5854A),  and  5  ^3/2  These  transitions  are  illustrated 

in  Fig.  21.  In  the  oast  the  common  procedure  has  been  to  observe  only  the 
strong  resonance  line,  A554A,  and  occasionally  the  other  member  of  the 
pair,  4934A,  to  obtain  column  densities  of  Ba'*‘.  A  radiative-transfer 
analysis  still  in  preparation  at  Los  Alamos  Scientific  Laboratory  by  Horak 
and  l^hitaker  (1980)  shows  that  observations  in  the  resonance  lines  alone 
can  give  erroneous  values  of  the  Integrated  density,  but  that  the 
observations  can  be  interpreted  unambiguously  if  supplemented  by 
simultaneous  photometry  of  the  fluorescent  lines. 

The  effect  is  -ne  commonly  known  as  the  curve  of  growth  and  is 
illustrated  in  Fig.  22.  At  column  densities  of  Ba+  of  N  *  10^^'  cm  the 
4554A  line  begins  to  become  saturated.  A  correction  is  not  simply  made  for 
the  loss  of  radiation  through  attenuation,  since  some  of  the  loss  is 
restored  through  secondary  scattering  back  into  the  line  of  sight.  The 
only  prnctical  means  of  ascertaining  the  column  density  is  through 
simultaneous  observations  of  the  lines  6497,  5854,  and/or  6142A,  which  have 
diminished  populations  in  their  lower  levels  and  consequently  are  slower  to 
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Figure  22  Curve  of  growth  for  one  resonance  transition  (4554A)  and  two  fluorescent  transitions 
(5854A  and  6142A).  Note  that  the  resonance  transition  saturates  at  high  column  den¬ 
sities  of  Ba**^  ions.  Thus  observations  at  4554A  cannot  be  used  to  correctly  measure 
Ba^  column  densities  above  ^^2  x  10^^  cm”^,  whereas  observations  at  the  other  two 
wavelengths  could  extend  measurements  by  at  least  an  order  of  magnitude.  After  Horak 
and  Whitaker  (1980). 
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saturate.  It  can  be  seen  in  Fig.  22  that  even  6142A  Is  beginning  to 
saturate  at  column  densities  of  10l2Ba+  cm“2  for  the  viewing  geometry 
indicated  on  the  figure.  So  long  as  the  cloud  remains  optically  thin  in 
the  resonance  lines,  there  is  no  ambiguity,  if  the  spectra  are  "pure." 

However,  another  problem  that  has  arisen  in  the  past  and  which 
further  argues  for  some  redundancy  in  the  spectral  observations  is  the 
leakage  of  filters  to  undeslred  emission.  Particularly,  the  4607A  line  of 
Sr  is  strong  and  seems  to  be  an  inevitable  contaminant  of  barium.  Unlike 
barium  the  strontium  is  essentially  unlonlzable,  and  the  4607A  line  may 
penetrate  the  4554A  filter. 

Because  the  forthcoming  "PLACES"  experiments  will  produce  optical 
thicknesses  in  the  range  40-60,  it  is  important  that  the  curve-of  growth 
effects  be  anticipated.  It  should  be  toted  that  curve-of-growth  effects 
depend  strongly  on  the  hyperfine  structure  of  the  scattering  atom  or  ion 
(J.  W.  Chamberlain,  1960).  In  the  LASL  work  of  Horak  and  Whitaker,  the 
complex  array  of  HFS  from  the  variety  of  Isotopes  has  been  accounted  for. 
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